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Abstract. The considerable part of the parameter space in the MSSM corresponding to 
the infrared quasi fixed point scenario is excluded by LEP II bounds on the lightest Higgs 
boson mass. In the NMSSM the mass of the lightest Higgs boson reaches its maximum 
value in the strong Yukawa coupling limit when Yukawa couplings are essentially larger 
than gauge ones at the Grand Unification scale. In this case the renormalization group 
flow of Yukawa couplings and soft SUSY breaking terms is investigated. The quasi-fixed 
and invariant lines and surfaces are briefly discussed. The coordinates of the quasi-fixed 
points, where all solutions are concentrated, are given. 

1 Introduction 

The search for the Higgs boson remains one of the top priorities for existing ac- 
celerators as well as for those still at the design stage. This is because this boson 
plays a key role in the Standard Model which describes all currently available 
experimental data with a high degree of accuracy. As a result of the spontaneous 
symmetry breaking SU(2) <g> U(1) the Higgs scalar acquires a nonzero vacuum 
expectation value without destroying the Lorentz invariance, and generates the 
masses of all fermions and vector bosons. An analysis of the experimental data 
using the Standard Model has shown that there is a 95% probability that its mass 
will not exceed 210 GeV |1 [. At the same time, assuming that there are no new 
fields and interactions and also no Landau pole in the solution of the renormal- 
ization group equations for the self-action constant of Higgs fields up to the scale 
Mpi « 2.4 • 1 1 8 GeV, we can show that m h < 1 80 GeV 0,O1- In this case, phys- 
ical vacuum is only stable provided that the mass of the Higgs boson is greater 
than 135 GeV |2[-|6[. However, it should be noted that this simplified model does 
not lead to unification of the gauge constants 1 7 1 and a solution of the hierarchy 
problem |8[. As a result, the construction of a realistic theory which combines all 
the fields and interactions is extremely difficult in this case. 

Unification of the gauge constants occurs naturally on the scale Mx ~ 3 • 
10 16 GeV within the supersymmetric generalisation of the Standard Model, i.e., 
the Minimal Supersymmetric Standard Model (MSSM) |7|. In order that all the 

* Editors' note: This contribution was intended for the Holger Bech Nielsen's Festschrift 
(Vol. 1 of this Proceedings), but was received late, so we include it here. 
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fundamental fermions acquire mass in the MSSM, not one but two Higgs dou- 
blets Hi and H2 must be introduced in the theory, each acquiring the nonzero 
vacuum expectation value V] and V2 where v 2 = v 2 + v\ — (246 GeV) 2 . The spec- 
trum of the Higgs sector of the MSSM contains four massive states: two CP-even, 
one CP-odd, and one charged. An important distinguishing feature of the super- 
symmetric model is the existing of a light Higgs boson in the CP-even sector. 
The upper bound on its mass is determined to a considerable extent by the value 
tan |3 = V2/V1 . In the tree-level approximation the mass of the lightest Higgs bo- 
son in the MSSM does not exceed the mass of the Z-boson (Mz ~ 91.2 GeV): 
m-h < M-zl cos 2 |3 |9|. Allowance for the contribution of loop corrections to the 
effective interaction potential of the Higgs fields from a t-quark and its super- 
partners significantly raises the upper bound on its mass: 



m H < a/M 2 cos 2 2(3 + A. 



(1) 



Here A are the loop corrections |10|,| 11 [. The values of these corrections are pro- 
portional to where rat is the running mass of t-quark which depends log- 
arithmically on the supersymmetry breaking scale Ms and is almost indepen- 
dent of the choice of tan (3. In |3[,|5[,|6| bounds on the mass of the Higgs boson 
were compared in the Minimal Standard and Supersymmetric models. The up- 
per bound on the mass of the light CP-even Higgs boson in the MSSM increases 
with increasing tan (3 and for tan |3 ^> 1 in realistic supersymmetric models with 
M s < 1 000 GeV reaches 1 25 - 1 28 GeV. 

However, a considerable fraction of the solutions of the system of MSSM 
renormalization group equations is focused near the infrared quasi-fixed point 
at tan |3 ~ 1 . In the region of parameter space of interest to us (tan (3 <C 50) the 
Yukawa constants of a b-quark (lib) and a T-lepton (h T ) are negligible so that an 
exact analytic solution can be obtained for the one-loop renormalization group 
equations [12J. For the Yukawa constants of a t-quark H t (t) and the gauge con- 
stants gi(t) its solution has the following form: 
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where the index i has values between 1 and 3, 
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The variable t is determined by a standard method t = ln( M^ /q 2 ) . The boundary 
conditions for the renormalization group equations are usually set at the grand 
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unification scale Mx (t = 0) where the values of all three Yukawa constants are 
the same: &i (0) = a.z{0) — ocs(0) — ao. On the electroweak scale where h^(0) ^> 1 
the second term in the denominator of the expression describing the evolution of 
Yt(t) is much smaller than unity and all the solutions are concentrated in a nar- 
row interval near the quasi-fixed point YQpp(t) = E(t)/6F(t) 1131 . In other words 
in the low-energy range the dependence of Y t (t) on the initial conditions on the 
scale Mx disappears. In addition to the Yukawa constant of the t-quark, the cor- 
responding trilinear interaction constant of the scalar fields A t and the combina- 
tion of the scalar masses Til — rag + + m| also cease to depend on A t (0) and 
9Jtt (0) as Y t (0) increases. Then on the electroweak scale near the infrared quasi- 
fixed point A t (t) and QJt^ (t) are only expressed in terms of the gaugino mass 
on the Grand Unification scale. Formally this type of solution can be obtained if 
Y t (0) is made to go to infinity. Deviations from this solution are determined by 
ratio l/6F(t)Y t (0) which is of the order of 1 /10h^(0) on the electroweak scale. 

The properties of the solutions of the system of MSSM renormalization group 
equations and also the particle spectrum near the infrared quasi-fixed point for 
tan |3 ~ 1 have been studied by many authors 1 14 1,| 15 1. Recent investigations 1 15 1- 
|17| have shown that for solutions Y t (t) corresponding to the quasi-fixed point 
regime the value of tan |3 is between 1 .3 and 1 .8. These comparatively low val- 
ues of tan |3 yield significantly more stringent bounds on the mass of the light- 
est Higgs boson. The weak dependence of the soft supersymmetry breaking pa- 
rameters At(t) and (t) on the boundary conditions near the quasi-fixed point 
means that the upper bound on its mass can be calculated fairly accurately. A the- 
oretical analysis made in |15|,|16| showed that ran does not exceed 94 ± 5 GeV. 
This bound is 25 — 30 GeV below the absolute upper bound in the Minimal Super- 
symmetric Model. Since the lower bound on the mass of the Higgs boson from 
LEP II data is 1 1 3 GeV 1 1 |, which for the spectrum of heavy supersymmetric par- 
ticles is the same as the corresponding bound on the mass of the Higgs boson in 
the Standard Model, a considerable fraction of the solutions which come out to 
a quasi-fixed point in the MSSM, are almost eliminated by existing experimental 
data. This provides the stimulus for theoretical analyses of the Higgs sector in 
more complex supersymmetric models. 

The simplest expansion of the MSSM which can conserve the unification of 
the gauge constants and raise the upper bound on the mass of the lightest Higgs 
boson is the Next-to-Minimal Supersymmetric Standard Model (NMSSM) |18[- 
|20|. By definition the superpotential of the NMSSM is invariant with respect to 
the discrete transformations y' a — e 27I1 / 3 u a of the Z3 group |19[ which means 
that we can avoid the problem of the |j.-term in supergravity models. The thing is 
that the fundamental parameter a should be of the order of M P i since this scale is 
the only dimensional parameter characterising the hidden (gravity) sector of the 
theory. In this case, however, the Higgs bosons Hi and H2 acquire an enormous 
mass Hl ~~ \± 2 ~~ Mj^ and no breaking of SU(2) ® U(1) symmetry occurs. 
In the NMSSM the term (o.(Ai A2) in the superpotential is not invariant with re- 
spect to discrete transformations of the Z3 group and for this reason should be 
eliminated from the analysis = 0). As a result of the multiplicative nature of 
the renormalization of this parameter, the term u.(q) remains zero on any scale 
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q < Mx -j- Mpi. However, the absence of mixing of the Higgs doublets on elec- 
tro weak scale has the result that H i acquires no vacuum expectation value as a 
result of the spontaneous symmetry breaking and d-type quarks and charged 
leptons remain massless. In order to ensure that all quarks and charged leptons 
acquire nonzero masses, an additional singlet superfield Y with respect to gauge 
SU(2) (gi U(1 ) transformations is introduced in the NMSSM. The superpotential 
of the Higgs sector of the Nonminimal Supersymmetric Model |18|-|20| has the 
following form: 

Wh = AY(Ai A2) + y Y 3 . 0) 

As a result of the spontaneous breaking of SU(2) ® U(1) symmetry, the field Y 
acquires a vacuum expectation value ((Y) = y/>/2) and the effective u-term (u. 
Ay / y/2) is generated. 

In addition to the Yukawa constants A and x, and also the Standard Model 
constants, the Nonminimal Supersymmetric Model contains a large number of 
unknown parameters. These are the so-called soft supersymmetry breaking pa- 
rameters which are required to obtain an acceptable spectrum of superpartners of 
observable particles form the phenomenological point of view. The hypothesis on 
the universal nature of these constants on the Grand Unification scale allows us 
to reduce their number in the NMSSM to three: the mass of all the scalar particles 
mo, the gaugino mass and the trilinear interaction constant of the scalar 

fields A. In order to avoid strong CP-violation and also spontaneous breaking of 
gauge symmetry at high energies (M P i 2> E ^> m t ) as a result of which the scalar 
superpartners of leptons and quarks would require nonzero vacuum expectation 
values, the complex phases of the soft supersymmetry breaking parameters are 
assumed to be zero and only positive values of tRq are considered. Naturally uni- 
versal supersymmetry breaking parameters appear in the minimal supergravity 
model |23| and also in various string models |22[,|24[. In the low-energy region 
the hypothesis of universal fundamental parameters allows to avoid the appear- 
ance of neutral currents with flavour changes and can simplify the analysis of 
the particle spectrum as far as possible. The fundamental parameters thus de- 
termined on the Grand Unification scale should be considered as boundary con- 
ditions for the system of renormalization group equations which describes the 
evolution of these constants as far as the electroweak scale or the supersymmetry 
breaking scale. The complete system of the renormalization group equations of 
the Nonminimal Supersymmetric Model can be found in |25|, |26|. These exper- 
imental data impose various constraints on the NMSSM parameter space which 
were analysed in |27|,|28|. 

The introduction of the neutral field Y in the NMSSM potential leads to the 
appearance of a corresponding F-term in the interaction potential of the Higgs 
fields. As a consequence, the upper bound on the mass of the lightest Higgs boson 
is increased: 

m H < yyV 2 sin 2 2|3 + Micos 2 2|3 + A<V + a\] ] . (4) 

The relationship was obtained in the tree-level approximation (An = 0) in 
|20|. However, loop corrections to the effective interaction potential of the Higgs 
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fields from the t-quark and its superpartners play a very significant role. In 
terms of absolute value their contribution to the upper bound on the mass of 
the Higgs boson remains approximately the same as in the Minimal Supersym- 
metric Model. When calculating the corrections A-j^ and A^ we need to re- 
place the parameter (a. by Xy/ y/2. Studies of the Higgs sector in the Nonmini- 
mal Supersymmetric model and the one-loop corrections to it were reported in 
1211,1251,1281-1311. In |6| the upper bound on the mass of the lightest Higgs bo- 
son in the NMSSM was compared with the corresponding bounds on ran in the 
Minimal Standard and Supersymmetric Models. The possibility of a spontaneous 
CP-violation in the Higgs sector of the NMSSM was studied in I5TI.153. 

It follows from condition J3J that the upper bound on mn increases as A in- 
creases. Moreover, it only differs substantially from the corresponding bound in 
the MSSM in the range of small tan (3 . For high values (tan P » 1) the value 
of sin2|3 tends to zero and the upper bounds on the mass of the lightest Higgs 
boson in the MSSM and NMSSM are almost the same. The case of small tan |3 is 
only achieved for fairly high values of the Yukawa constant of a t-quark h t on the 
electroweak scale (h t (to) > 1 where to = lnfM^/m^ )), and tan (3 decreases with 
increasing h t (to). However, an analysis of the renormalization group equations 
in the NMSSM shows that an increase of the Yukawa constants on the electroweak 
scale is accompanied by an increase of h t (0) and A(0) on the Grand Unification 
scale. It thus becomes obvious that the upper bound on the mass of the lightest 
Higgs boson in the Nonminimal Supersymmetric model reaches its maximum on 
the strong Yukawa coupling limit, i.e., when h t (0) > gi(0) and A(0) > gi(0). 

2 Renormalization of the Yukawa couplings 

From the point of view of a renormalization group analysis, investigation of the 
NMSSM presents a much more complicated problem than investigation of the 
minimal SUSY model. The full set of renormalization group equations within 
the NMSSM can be found in |25|,|26|. Even in the one-loop approximation, this 
set of equations is nonlinear and its analytic solution does not exist. All equa- 
tions forming this set can be partitioned into two groups. The first one con- 
tains equations that describe the evolution of gauge and Yukawa coupling con- 
stants, while the second one includes equations for the parameters of a soft break- 
down of SUSY, which are necessary for obtaining a phenomenologically accept- 
able spectrum of superpartners of observable particles. Since boundary condi- 
tions for three Yukawa coupling constants are unknown, it is very difficult to per- 
form a numerical analysis of the equations belonging to the first group and of the 
full set of the equations given above. In the regime of strong Yukawa coupling, 
however, solutions to the renormalization group equations are concentrated in a 
narrow region of the parameter space near the electroweak scale, and this consid- 
erably simplifies the analysis of the set of equations being considered. 

In analysing the nonlinear differential equations entering into the first group, 
it is convenient to use the quantities p t , Pa, p*, Pi , and P2, defined as follows: 

M Y t (t) ... Y A (t) Y K {t) 

Pt t )=—. r- , P\[t) = —-, P*(t)=— : — , 

a 3 (t) a 3 (t] a 3 (t] 
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ait 
Pi (t) = — 



P2(t) 



where oci(t) 

^ 2 (t)/(47t) 2 . 
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Fig. 1. The values of the Yukawa couplings at the electroweak scale corresponding to the 
initial values at the GUT scale uniformly distributed in a square 2 < h. 2 (0) , A 2 (0) < 10. The 
thick and thin curves represent, respectively, the invariant and the Hill line. The dashed 
line is a fit of the values (p t (to), Pa (to)) for 20 < h 2 (0),A 2 (0) < 100. 



Let us first consider the simplest case of x — 0. The growth of the Yukawa 
coupling constant A(t ) at a fixed value of h t (t ) results in that the Landau pole 
in solutions to the renormalization group equations approaches the Grand Uni- 
fication scale from above. At a specific value A (to) = A max , perturbation theory 
at q - Mx cease to be applicable. With increasing (decreasing) Yukawa coupling 
constant for the b-quark, A max decreases (increases). In the (pt, Pa) plane, the de- 
pendence A^^fh 2 ) is represented by a curve bounding the region of admissible 
values of the parameters Pt(to) and Pa (to). At Pa = 0, this curve intersects the 
abscissa at the point pt = p^ FP (to). This is the way in which there arises, in the 
( Pt , Pa ) plane, the quasi-fixed (or Hill) line near which solutions to the renormal- 
ization group equations are grouped (see Fig. 1). With increasing A 2 (0) and h 2 (0), 
the region where the solutions in questions are concentrated sharply shrinks, and 
for rather large initial values of the Yukawa coupling constants they are grouped 
in a narrow stripe near the straight line 



pt(to) +0.506p A (t ) = 0.91, 



(5) 
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Fig. 2. The values of the Yukawa couplings at the electroweak scale corresponding to the 
initial values at the GUT scale uniformly distributed in a square 20 < h 2 (0) , A 2 (0) < 1 00. 
The dashed line is a fit of the values (pt(to), Pa (to)) for 20 < h 2 (0), A 2 (0) < 100. 



which can be obtained by fitting the results of numerical calculations (these re- 
sults are presented in Fig. 2). Moreover, the combination h 2 (to) + 0.506A 2 (to) of 
the Yukawa coupling constants depends much more weakly on A 2 (0) and h 2 (0) 
than A 2 (to ) and h 2 (to ) individually 1 36 [. In other words, a decrease in A 2 (to ) com- 
pensates for an increase in h 2 (to), and vice versa. The results in Fig. 3, which 
illustrate the evolution of the above combinations of the Yukawa coupling con- 
stants, also confirm that this combination is virtually independent of the initial 
conditions. 

In analysing the results of numerical calculations, our attention is engaged 
by a pronounced nonuniformity in the distribution of solutions to the renormal- 
ization group equations along the infrared quasi-fixed line. The main reason for 
this is that, in the regime of strong Yukawa coupling, the solutions in question are 
attracted not only to the quasi-fixed but also to the infrared fixed (or invariant) 
line. The latter connects two fixed points. Of these, one is an infrared fixed point 
of the set of renormalization group equations within the NMSSM (p t = 7/18, 
Pa = 0, pi = 0, p2 = 0) 1 33 1, while the other fixed point (px/pt — 1 ) corresponds 
to values of the Yukawa coupling constants in the region Y t , Yx ^> bi\, in which 
case the gauge coupling constants on the right-hand sides of the renormalization 
group equations can be disregarded |34[. For the asymptotic behaviour of the 
infrared fixed line at p t , p\ 3> 1 we have 



8 2 
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Fig. 3. Evolution of the combination p t (t) + 0.506pA (t) of the Yukawa couplings from the 
GUT scale (t = 0) to the electroweak scale (t = to ) for x 2 — and for various initial values 
h?(0)-Fig. 3a, A 2 (0) - Fig. 3b. 
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while in the vicinity of the point p t = 7/1 8, p A = we have 

p A ~(p t -7/18) 25 / U 

The infrared fixed line is invariant under renormalization group transforma- 
tions - that is, it is independent of the scale at which the boundary values Y t (0) 
and Y A (0) are specified and of the boundary values themselves. If the bound- 
ary conditions are such that Y t (0) and Y\(0) belong to the fixed line, the evolu- 
tion of the Yukawa coupling constants proceeds further along this line toward 
the infrared fixed point of the set of renormalization group equations within the 
NMSSM. With increasing t, all other solutions to the renormalization group equa- 
tions are attracted to the infrared fixed line and, for t/(4n) ^> 1, approach the 
stable infrared fixed point. From the data in Figs. 1 and 2, it follows that, with 
increasing Y t (0) and Y\(0), all solutions to the renormalization group equations 
are concentrated in the vicinity of the point of intersection of the infrared fixed 
and the quasi-fixed line: 

pf p (t )= 0.803, pf (to )= 0.224. 

Hence, this point can be considered as the quasi-fixed point of the set of renor- 
malization group equations within the NMSSM at k — 0. 

In a more complicated case where all three Yukawa coupling constants in 
the NMSSM are nonzero, analysis of the set of renormalization group equations 
presents a much more difficult problem. In particular, invariant (infrared fixed) 
and Hill surfaces come to the fore instead of the infrared fixed and quasi-fixed 
lines. For each fixed set of values of the coupling constants Y t (t ) and Y^(t ), 
an upper limit on Y\ (to ) can be obtained from the requirement that perturbation 
theory be applicable up to the Grand Unification scale Mx. A change in the val- 
ues of the Yukawa coupling constants fi t and x at the electroweak scale leads to 
a growth or a reduction of the upper limit on Y A (t ). The resulting surface in the 
(p t , p^, Pa) space is shown in Fig. 4. In the regime of strong Yukawa coupling, 
solutions to the renormalization group equations are concentrated near this sur- 
face. In just the same way as in the case of Y^ = 0, a specific linear combination 
of Y t , Y\, and Y^ is virtually independent of the initial conditions for Y|(0) — > oo: 

Pt(to) + 0.72p A (t ) + 0.33p„(t ) = 0.98. (6) 

The evolution of this combination of Yukawa couplings at various initial values 
of the Yukawa coupling constants is illustrated in Fig. 5. 

On the Hill surface, the region that is depicted in Fig. 4 and near which the 
solutions in question are grouped shrinks in one direction with increasing initial 
values of the Yukawa coupling constants, with the result that, at Y t (0), Y^(0), and 
Y A (0) ~ 1, all solutions are grouped around the line that appears as the result 
of intersection of the quasi-fixed surface and the infrared fixed surface, which 
includes the invariant lines lying in the p^ = and p A = planes and connect- 
ing the stable infrared point with, respectively, the fixed point p A /pt = 1 and 
the fixed point p^/Pt = 1 in the regime of strong Yukawa coupling. In the limit 
Pt , Pa , > 1 , in which case the gauge coupling constants can be disregarded, 
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Fig. 4. Quasi-fixed surface in the ( pt , p* , Pa ) space. The shaded part of the surface rep- 
resents the region near which the solutions corresponding to the initial values 2 < 
h. 2 t {0),x 2 {0),A 2 {0) < 10 -Fig. 4a, 20 < hl(0),x 2 (0),A 2 {0) < 1 00 - Fig. 4b are concen- 
trated. 
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Fig. 5. Evolution of the combination p t +0.720px +0.3330p^ of the Yukawa couplings from 
the GUT scale (t = 0) to the electroweak scale (t = to) for various initial values h|(0) - 
Fig. 5a, A 2 (0) - Fig. 5b, x 1 (0) - Fig. 5c. 
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the fixed points Px/pt = 1, P^/Pt = and p^Pt = 1, PA/pt = cease to be 
stable. Instead of them, the stable fixed point Ra = 3/4, R^ — 3/8 |34| appears 
in the (Ra, R^) plane, where Ra = PA/pt and R^ = p^/pt- In order to investigate 
the behaviour of the solutions to the renormalization group equations within the 
NMSSM, it is necessary to linearise the set of these equations in its vicinity and 
set ctt = 0. As a result, we obtain 




where Rao = Ra(0), R>,o = R*r(0), p t o = 

Pt(t) = - — ^ . From 171, it follows that the fixed point Ra 

I + 7ptot 

arises as the result of intersection of two fixed lines in the (Ra>R^) plane. The 

solutions are attracted most strongly to the line -Ra H pR>r = f 1 H — t= 

since A] ^ A2. This line passes through three fixed points in the (Ra, R>t) plane: 
(1,0), (3/4,3/8), and (0, 1). In the regime of strong Yukawa coupling, the fixed 
line that corresponds, in the ( Pt , P>o P a ) space, to the line mentioned immediately 
above is that which lies on the invariant surface containing a stable infrared fixed 
point. The line of intersection of the Hill and the invariant surface can be obtained 
by mapping this fixed line into the quasi-fixed surface with the aid of the set of 
renormalization group equations. For the boundary conditions, one must than 
use the values A 2 (0), x 2 (0), and h 2 (0) » 1 belonging to the aforementioned fixed 
line. 

In just the same way as infrared fixed lines, the infrared fixed surface is in- 
variant under renormalization group transformations. In the evolution process, 
solutions to the set of renormalization group equations within the NMSSM are 
attracted to this surface. If boundary conditions are specified n the fixed sur- 
face, the ensuing evolution of the coupling constants proceeds within this sur- 
face. To add further details, we not that, near the surface being studied and on it, 
the solutions are attracted to the invariant line connecting the stable fixed point 
(PA/Pt = 3/4, p>t/pt = 3/8) in the regime of strong Yukawa coupling with the 
stable infrared fixed point within the NMSSM. In the limit p t , p^, Pa ^> 1, the 
equation for this line has the form 

3 176 3 7 

P* = 4 Pt -4T7 + T39 P2 -4T7 Pl ' (H) 
3 56 18 68 ( ' 

P *-8 Pt ~4l7~T39 P2 ~ 2085 Pl - 
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As one approaches the infrared fixed point, the quantities Pa and tend to zero: 
Pa ~ (Pt — 7/18) 25/14 and p^ - (p t — 7/18) 9/7 . This line intersects the quasi-fixed 
surface at the point 

pf p (t ) = 0.82, p2 FP (t ) = 0.087, pf p (t ) = 0.178. 

Since all solutions are concentrated in the vicinity of this point for Y t (0), Y\(0), 
Yjt(0) — > oo, it should be considered as a quasi-fixed point for the set of renor- 
malization group equations within the NMSSM. We note, however, that the so- 
lutions are attracted to the invariant line and to the quasi-fixed line on the 
Hill surface. This conclusion can be drawn from the an analysis of the behaviour 
of the solutions near the fixed point (Ra = 3/4, R^ = 3/8) (see 0). Once the 

1„ 1 3/. 1 



solutions have approached the invariant line -Ra H — pR x = „ M H — t= )/ their 

evolution is governed by the expression (e(t)) 085 , where e(t) = p t (t)/p t o- This 
means that the solutions begin to be attracted to the quasi-fixed point and to the 
invariant line with a sizable strength only when Yi (0) reaches a value of 1 2 , at 
which perturbation theory is obviously inapplicable. Thus, it is not the infrared 
quasi-fixed point but the quasi-fixed line on the Hill surface (see Fig. 4) that, 
within the NMSSM, plays a key role in analysing the behaviour of the solutions 
to the renormalization group equations in the regime of strong Yukawa coupling, 
where all Yt(0) are much greater than ocq. 



3 Renormalization of the soft SUSY breaking parameters 

If the evolution of gauge and Yukawa coupling constants is known, the remaining 
subset of renormalization group equations within the MNSSM can be treated as 
a set of linear differential equations for the parameters of a soft breakdown of 
supersymmetry. For universal boundary conditions, a general solution for the 
trilinear coupling constants Ai(t) and for the masses of scalar fields rn 2 (t) has 
the form 

At(t) = et(t)A + f t (t)M 1/2 , (9) 
m 2 (t) = ot(t)Tng + bi(t)M 2 /2 + d(t)AM 1/2 + d t (t)A 2 . (10) 

The functions et(t), ft(t), a-t(t), bi(t), Ct(t), and di(t), which determine the evo- 
lution of A-i(t) and m 2 (t), remain unknown, since an analytic solution to the full 
set of renormalization group equations within the NMSSM is unavailable. These 
functions greatly depend on the choice of values for the Yukawa coupling con- 
stants at the Grand Unification scale Mx. At the electroweak scale t = to, rela- 
tions and J10I specify the parameters A? (to) and m 2 (to) of a soft breaking of 
supersymmetry as functions of their initial values at the Grand Unification scale. 
The results of our numerical analysis indicate that, with increasing Yj,(0), 

where Y t (t) = ^ ^] , Y A (t) = f * and Y*(t) = * , the functions ejto), 

Ci(to), and di(to) decrease and tend to zero in the limit Y;,(0) — > oo, relations 
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and (1101 becoming much simpler in this limit. Instead of the squares of the scalar 
particle masses, it is convenient to consider their linear combinations 

mt{ (t) = mj (t) + m 2 Q (t) + mu (t) , 

OnJ(t)=m?(t} + TnJ(t)+m*(t), (11) 
Tl 2 Jt)=3m 2 y (t) 

in analysing the set of renormalization group equations. In the case of universal 
boundary conditions, the solutions to the differential equations for 971? (t) can be 
represented in the same form as the solutions for m?(t) (see <10» ; that is 

St7t?(t) = 3ttiCt)-mg + 6 t (t)M? /2 + Ci(t)AM 1/2 + di(t)A 2 . (12) 

Since the homogeneous equations for AJt) and 971? (t) have the same form, the 
functions Qi(t) and et(t) coincide; in the limit of strong Yukawa coupling, the 
nig dependence disappears in the combinations 11 U of the scalar particle masses 
as the solutions to the renormalization group equations for the Yukawa coupling 
constants approach quasi-fixed points. This behaviour of the solutions implies 
that A-i(t) and 971? (t) corresponding to Yi(0) > cct(0) also approach quasi-fixed 
points. As we see in the previous section, two quasi-fixed points of the renor- 
malization group equations within the NMSSM are of greatest interest from the 
physical point of view. Of these, one corresponds to the boundary conditions 
Y t (0) = Y A (0) > &i(0) and Y^(0) = for the Yukawa coupling constants. The 
fixed points calculated for the parameters of a soft breaking of supersymmetry 
by using these values of the Yukawa coupling constants are 

pQ FP (t )«1.77, pQ F f(to) « 6.09, 

pQf (to) « -0.42, pQ FP (t ) a -2.28, ' 

where PA t (t) = Ai(t)/Mi/2 and p OT 2(t) = Wl 2 /M 2 /2 . Since the coupling con- 
stant h for the self-interaction of neutral scalar fields is small in the case being 
considered, A M (t) and 971 2 (t) do not approach the quasi-fixed point. Nonethe- 
less, the spectrum of SUSY particles is virtually independent of the trilinear cou- 
pling constant A^ since k — > 0. 

In just the same way, one can determine the position of the other quasi-fixed 
point for Ai(t) and 971? (t), that which corresponds to Rao = 3/4, R^o = 3/8. The 
results are 

pQ FP (t ) « 1.73, pQf(to) « -0.43, pQ FP (t ) w 0.033, 

/v, ^ /v« 

pQ F r(t )«6.02, pQ F 2 p (t ) « -2.34, pQ FP (t )«0.29, 



where R A0 = Y A (0)/Y t (0) and R^ = Y^(0)/Y t (0). It should be noted that, in the 
vicinities of quasi-fixed points, we have pSff (to) < 0. Negative values of 971? (to) 

lead to a negative value of the parameter m|(to) in the potential of interaction 
of Higgs fields. In other words, an elegant mechanism that is responsible for a 
radiative violation of SU(2) <£> U(1 ) symmetry and which does not require intro- 
ducing tachyons in the spectrum of the theory from the outset survives in the 
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regime of strong Yukawa coupling within the NMSSM. This mechanism of gauge 
symmetry breaking was first discussed in 1 35 [ by considering the example of the 
minimal SUSY model. 

By using the fact that 9Jt? (t) as determined for the case of universal boundary 
conditions is virtually independent of m.Q, we can predict at (to) values near the 
quasi-fixed points (see |37[). The results are 

1 ) Rao = 1 , Rxo = 0, 

1 4 5 

a y [t ) = Q u (to) = =, a! (t ) = a q (t ) = -, a 2 (t ) = -- ; 

' ' 1 (15) 

2) R A0 = 3/4, R^ = 3/8, 

2 5 1 

ay (to) = 0, ai (t ) = -a 2 (t ) = -, a q (t ) = -, a u (t ) = - . 

To do this, it was necessary to consider specific combinations of the scalar particle 
masses, such as — 2tuq , iUg + — m.^ + , and — 2mf (at k = 0), that 
are not renormalized by Yukawa interactions. As a result, the dependence of the 
above combinations of the scalar particle masses on at the electroweak scale 
is identical to that at the Grand Unification scale. The predictions in J15> agree 
fairly well with the results of numerical calculations. 



-0.3 



-0.4 



-0.5 



1.72 



1.74 



1.76 



1.78 



A t /M 1/: 



Fig. 6. The values of the trilinear couplings A t and Aa at the electroweak scale correspond- 
ing to the initial values uniformly distributed in the ( A t , Aa ) plane, calculated at x 1 — 
and H? (0) = A 2 (0) = 20. The straight line is a fit of the values (A t (t ), A A (t )). 



Let us now consider the case of nonuniversal boundary conditions for the 
soft SUSY breaking parameters. The results of our numerical analysis, which are 
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6.02 6.04 6.06 6.08 6.1 6.12 dJl^/My 2 

Fig. 7. The values of the combinations of masses DJl 2 and 9Jt A at the electroweak scale 
corresponding to the initial values uniformly distributed in the (dR^/M. 2 /2 ,^ytx/M. 2 /2 ) 
plane, calculated at = 0, h^(0) = A 2 (0) = 20, and A t (0) = A A (0) = 0. The straight line 
is a fit of the values [Til (to ) , Wll (t ) ). 



illustrated in Figs. 6 and 7, indicate that, in the vicinity of the infrared fixed point 
at = 0, solutions to the renormalization group equations at the electroweak 
scale are concentrated near some straight lines for the case where the simula- 
tion was performed by using boundary conditions uniformly distributed in the 
(A t , A A ) and the (501^,971^) plane. The strength with which these solutions are 
attracted to them grows with increasing Yi(0). The equations for the lines being 
considered can be obtained by fitting the numerical results displayed in Figs. 6 
and 7. This yields 

A t + 0.147A A = 1.70Mi /2 , 

M\ + 0.14750^ = 5.76Mf /2 . ^ 

For Yjt(0) 3> oco solutions to the renormalization group equations are grouped 
near planes in the space of the parameters of a soft breaking of supersymmetry 
(A t , A X ,A X ) and (2K?,0rt£, Mi ) (see Figs. 8-10): 

A t +0.128A A +0.022A„ = 1.68M 1/2 , 

(17) 

ml + 0.1289Jt A + 0.022371* = 5.77M^ /2 . 

It can be seen from Figs. 8 and 9 that, as the values of the Yukawa coupling con- 
stants at the Grand Unification scale are increased, the areas of the surfaces near 
which the solutions At(t) and SPl?(t) are concentrated shrink in one of the direc- 
tions, with the result that, at Yt(0) ~ 1, the solutions to the renormalization group 
equations are attracted to one of the straight lines belonging to these surfaces. 
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Fig.8. Planes in the parameter spaces (At/M.1/2, A^/M.i/2, A^/M.1/2) - Fig. 8a, and 
(9J^/Mj /2 ,97^/Mf /2 ,9Jti/Mf /2 ) - Fig. 8b. The shaded parts of the planes correspond 
to the regions near which the solutions at tit (0) = 16, A 2 (0) = 12, and x 2 (0) = 6 are 
concentrated. The initial values At (0) and 9Jt 2 (0) vary in the ranges — M] / 2 < A < /2 
and < OT 2 (0) < 3M 2 /2 , respectively. 
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Fig. 9. Planes in the parameter spaces (At /Mi / 2 , Aa/M, /2) A^/M, /2 ) - Fig. 9a, and 
(9Jt?/M? /2) £Dt£/Mf /2l 9j£/M? /2 ) - Fig. 9b. The shaded parts of the planes correspond 
to the regions near which the solutions at hi (0) = 32, A 2 (0) = 24, and x 2 (0) = 12 are 
concentrated. The initial values At (0) and SOT? (0) vary in the ranges — M] /2 < A < M.^ / 2 
and < m\ (0) < 3M^ /2 , respectively. 



124 R. B. Nevzorov, K. A. Ter-Martirosyan and M. A. Trusov 



The numerical calculations also showed that, with increasing Yt(0), only in 
the regime of infrared quasi-fixed points (that is, at Rao = 1 , R>rO = or at Rao = 
3/4, R^o = 3/8) ei(to) and at (to) decrease quite fast, in proportion to 1/Yt(0). 
Otherwise, the dependence on A and m.Q disappears much more slowly with 
increasing values of the Yukawa coupling constants at the Grand Unification scale 
-specifically, in proportion to (Yt(0))~ 6 , where 6 < 1 (for example, 6 =0.35 — 0.40 
at x = 0). In the case of nonuniversal boundary conditions, only when solutions 
to the renormalization group equations approach quasi-fixed points are these 
solutions attracted to the fixed lines and surfaces in the space of the parameters 
of a soft breaking of supersymmetry, and in the limit Yt(0) — > oo, the parameters 
A-t(t) and 971? (t) cease to be dependent on the boundary conditions. 

For the solutions of the renormalization group equations for the soft SUSY 
breaking parameters near the electroweak scale in the strong Yukawa coupling 
regime one can construct an expansion in powers of the small parameter e t (t) = 
V«(t)/Y,(0): 

;^ ( , t ) ) )=L^v ll (.)(^)(e l (t))- + ..., ,18, 

where at and |3t are constants of integration that can be expressed in terms of 
At(0) and 971? (0). The functions Vy (t) are weakly dependent on the Yukawa cou- 
pling constants at the scale Mx, and (0) = 1 . They appear upon renormalizing 
the parameters of a soft breaking of supersymmetry from q ~ 1 1 2 — 1 1 3 GeV 
to q - m.t. In equations <18> , we have omitted terms proportional to ?VL 1 /2, M.^ , 2/ 
Ai(0)Mi /2/ andAi(0)A,(0). 

At k — 0, we have two eigenvalues and two corresponding eigenvectors: 

3/ 7 ) ' U= (!-3 

whose components specify ( A t , A A ) and (9Jtf , 971^ ) ■ With increasing Y t (0) ~ Y A (0), 
the dependence on oco and |3o becomes weaker and the solutions at t = to are 
concentrated near the straight lines (At (<X1 ) , Aa (<xi ) ) and (97?^ (Pi), 971^ 1 ) )• In 
order to obtain the equations for these straight lines, it is necessary to set Aa(0) = 
-3A t (0) and 971^(0) = -3971^(0) at the Grand Unification scale. At the elec- 
troweak scale, there then arise a relation between At (to) and Aa (to ) and a relation 
between 97^ (to) and 97t^(t ): 

A t + 0.137A A = 1.70MV2, 

(19) 

97^ + 0.13797^ = 5.76M? /2 . 

These relations agree well with the equations deduced for the straight lines at 
Yt(0) - 1 by fitting the results of the numerical calculations | |16> . 

When the Yukawa coupling constant n is nonzero, we have three eigenvalues 
and three corresponding eigenvectors: 





I 1 \ 




A = 


3 + V5 | 

9 


■-( 
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-0.12 -0.1 -0.08 -0.06 -0.04 -0.02 

0.1278(A A /M 1/2 ) + 0.0223(A„/Mi /2 ) 

(a) 




-0.4 -0.375 -0.35 -0.325 -0.3 -0.275 -0.25 -0.225 



0.1278(SOTj/M 1 2 /2 ) + 0.0223(97^/71^) 

(b) 

Fig. 10. Set of points in planes (0.0223(A„/M.i /2 ) + 0.1278(Aa/M V2 ), A t /M 1/2 ) -Fig. 
10a, and (0.0223(5W 2 /M 2 /2 ) + 0.1278(OT£/M 2 /2 ), <H 2 /M 2 /2 ] - Fig. 10b, corresponding 
to the values of parameters of soft SUSY breaking for h? (0) = 32, A 2 (0) = 24, ;* 2 (0) = 
12, and for a uniform distribution of the boundary conditions in the parameter spaces 
(At,Ax,A„) and [ml,Ml,9Jti). The initial values A t (0) and 5W 2 (0) vary in the ranges 
-M,/2 <A<M, /2 and < (0) < 3M 2 /2 , respectively. The straight lines in Figs. 10a 
and 10b correspond to the planes in Figs. 9a and 9b, respectively. 
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whose components specify (A t , A A , A^) and (97tt, 971^,971^). An increase in Y A (0) ~ 
3 

2Y^(0) ~ ^Yt(O) leads to the following: first, the dependence of Ai(t] and 971? (t) 

on oco and (3o disappears, which leads to the emergence of planes in the space 
spanned by the parameters of a soft breaking of supersymmetry: 

A t + 0.103A A +0.0124A„ = 1.69M 1/2 , 
Ml + 0.10397t A + 0.0124971 2 , = 5.78M? /2 . (-0) 

After that, the dependence on cti and |3 1 becomes weaker at Yt (0) ~ 1 . This means 
that, with increasing initial values of the Yukawa coupling constants, solutions to 
the renormalization group equations are grouped near some straight lines and we 
can indeed see precisely this pattern in Figs. 8-10. All equations presented here 
for the straight lines and planes in the 971? space were obtained at A-t(O) =0. 

From relations fi"9l and l2"0l . it follows that At (to) and 971^ (to) are virtually 
independent of the initial conditions; that is, the straight lines and planes are 
orthogonal to the A t and 971^ axes. On the other hand, the A^(to) and 971^ (to) 
values that correspond to the Yukawa self-interaction constant Y^ for the neutral 
fields are fully determined by the boundary conditions for the parameters of a 
soft breaking of supersymmetry. For this reason, the planes in the (At, A A , A x ) 
and (97tt , 97t A , 97t?J spaces are virtually parallel to the A^ and 971^ axes. 

4 Conclusions 

In the strong Yukawa coupling regime in the NMSSM, solutions to the renormal- 
ization group equations for Y^(t) are attracted to quasi-fixed lines and surfaces 
in the space of Yukawa coupling constants and specific combinations of Pi(t) 
are virtually independent of their initial values at the Grand Unification scale. 
For Yi(0) — > oo, all solutions to the renormalization group equations are concen- 
trated near quasi-fixed points. These points emerge as the result of intersection 
of Hill lines or surfaces with the invariant line that connects the stable fixed point 
for Yi 3> oti with the stable infrared fixed point. For the renormalization group 
equations within the NMSSM, we have listed all the most important invariant 
lines and surfaces and studied their asymptotic behaviour for Y^ 3> at and in the 
vicinity of the infrared fixed point. 

With increasing Yj,(0), the solutions in question approach quasi-fixed points 
quite slowly; that is, the deviation is proportional to (et(t)) 6 , where e t (t) = 
Y t (t)/Y t (0) and 6 is calculated by analysing the set of the renormalization group 
equations in the regime of strong Yukawa coupling. As a rule, 5 is positive and 
much less than unity. By way of example, we indicate that, in the case where all 
three Yukawa coupling constants differ from zero, 6 0.085. Of greatest impor- 
tance in analysing the behaviour of solutions to the renormalization group equa- 
tions within the NMSSM at Y t (0), Y A (0), Y^(0) ~ 1 is therefore not the infrared 
quasi-fixed point but the line lying on the Hill surface and emerging as the in- 
tersection of the Hill and invariant surface. This line can be obtained by mapping 
the fixed points (1,0), (3/4, 3/8), and (0, 1 ) in the (R A , plane for Y t > cti into 
the quasi-fixed surface by means of renormalization group equations. 
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While Yt(t) approach quasi-fixed points, the corresponding solutions for the 
trilinear coupling constants AJt) characterising scalar fields and for the combi- 
nations 971? (t) of the scalar particle masses (see fTi) l cease to be dependent on 
their initial values at the scale Mx and, in the limit Yt(O) — > oo, also approach the 
fixed points in the space spanned by the parameters of a soft breaking of super- 
symmetry. Since the set of differential equations for AJt) and m?(t) is linear, the 
A, Mi /i, and ttlq dependence of the parameters of a soft breaking of supersym- 
metry at the electroweak scale can be explicitly obtained for universal boundary 
conditions. It turns out that, near the quasi-fixed points, all At(t) and all 971? (t) 
are proportional to /i and My 2 ' respectively Thus, we have shown that, in 
the parameter space region considered here, the solutions to the renormalization 
group equations for the trilinear coupling constants and for some combinations 
of the scalar particle masses are focused in a narrow interval within the infrared 
region. Since the neutral scalar field Y is not renormalized by gauge interactions, 
A x (t) and 971^ (t) are concentrated near zero; therefore they are still dependent 
on the initial conditions. The parameters At (to) and 971^ (to) show the weakest 
dependence on A and m.^ because these parameters are renormalized by strong 
interactions. By considering that the quantities 971? (to) are virtually independent 
of the boundary conditions, we have calculated, near the quasi-fixed points, the 
values of the scalar particle masses at the electroweak scale. 

In the general case of nonuniversal boundary conditions, the solutions to 
the renormalization group equations within the NMSSM for A-Jt) and 971? (t) are 
grouped near some straight lines and planes in the space spanned by the param- 
eters of a soft breaking of supersymmetry. Moving along these lines and surfaces 
as Yi(0) is increased, the trilinear coupling constants and the above combinations 
of the scalar particle masses approach quasi-fixed points. However, the depen- 
dence of these couplings on At(0) and 971? (0) dies out quite slowly, in proportion 
to (e t (t)) A , where A is a small positive number; as a rule, A < 1, For example, 
A = 3/7 at Y^ = and A w 0.0085 at Y^ ^ 0. The above is invalid only for the 
solutions Ai(t) and 971? (t) that correspond to universal boundary conditions for 
the parameters of a soft breaking of supersymmetry and to the initial values of 
Rao = 1 , Rxo = and Rao = 3/4, R^o = 3/8 for the Yukawa coupling con- 
stants at the Grand Unification scale. They correspond to quasi-fixed points of 
the renormalization group equations for Yt(t). As the Yukawa coupling constants 
are increased, such solutions are attracted to infrared quasi-fixed points in pro- 
portion to e t (t). 

Straight lines in the (A t , A A , A*} and (971^,971^,971^) spaces play a key role in 
the analysis of the behaviour of solutions for A-t(t) and 971? (t) in the case where 
Y t (0), Ya(0), Y^O) ~ 1. In the space spanned by the parameters of a soft break- 
ing of supersymmetry, these straight lines lie in the planes near which AJt) and 
971? (t) are grouped in the regime of strong Yukawa coupling at the electroweak 
scale. The straight lines and planes that were obtained by fitting the results of nu- 
merical calculations are nearly orthogonal to the A t and 97?^ axes. This is because 
the constants At (to) and 97tf (to) are virtually independent of the initial condi- 
tions at the scale Mx - On the other hand, the parameters A^(to) and 971^ (to) are 
determined, to a considerable extent, by the boundary conditions at the scale Mx- 
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At R A0 = 3/4 and R^o = 3/8, the planes in the (A t , A*, A*) and (SD^.SDt^SDt^) 
spaces are therefore parallel to the A^ and DJt^ axes. 
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The simplest effective dynamics describing the confinement mechanism in 
the pure gauge lattice U(l) theory is the dual Abelian Higgs model of scalar 
monopoles [1-3]. 

In the previous papers [4-6] the calculations of the U(l) phase transition (crit- 
ical) coupling constant were connected with the existence of artifact monopoles 
in the lattice gauge theory and also in the Wilson loop action model 1 6 1. In Ref .| 6 1 
we (L.V.L. and H.B.N.) have put forward the speculations of Refs.[4,5] suggesting 
that the modifications of the form of the lattice action might not change too much 
the phase transition value of the effective continuum coupling constant. In 1 6 1 the 
Wilson loop action was considered in the approximation of circular loops of radii 
R > a. It was shown that the phase transition coupling constant is indeed approx- 
imately independent of the regularization method: a cr it ~ 0.204, in correspon- 
dence with the Monte Carlo simulation result on lattice |7|: a cr tt ~ 0.20 ± 0.015. 

But in Ref s. [2,3] instead of using the lattice or Wilson loop cut-off we have 
considered the Higgs Monopole Model (HMM) approximating the lattice artifact 
monopoles as fundamental pointlike particles described by the Higgs scalar field. 

1 The Coleman-Weinberg effective potential for the Higgs 
monopole model 

The dual Abelian Higgs model of scalar monopoles (shortly HMM), describing 
the dynamics of confinement in lattice theories, was first suggested in Ref.| 1 1, and 
considers the following Lagrangian: 

L = -^F^(B) + l|(3^-iB li )0| 2 -U(0) ) where U(<D) = ^ 2 |0| 2 + ^|d>| 4 

(1) 

is the Higgs potential of scalar monopoles with magnetic charge g, and B u is the 
dual gauge (photon) field interacting with the scalar monopole field O. In this 
model A is the self-interaction constant of scalar fields, and the mass parameter 
ll 2 is negative. In Eq.Q the complex scalar field O contains the Higgs (4>) and 
Goldstone (x) boson fields: 

<& = * + iX- (2) 



132 



L.V. Laperashvili, D.A. Ryzhikh and H.B. Nielsen 



The effective potential in the Higgs Scalar ElectroDynamics (HSED) was first cal- 
culated by Coleman and Weinberg 1 8 1 in the one-loop approximation. The gen- 
eral method of its calculation is given in the review |9[. Using this method, we 
can construct the effective potential for HMM. In this case the total field system 
of the gauge (B,^) and magnetically charged (O) fields is described by the parti- 
tion function which has the following form in Euclidean space: 



Z = 



[DB][DO][DO+] e 



-s 



(3) 



where the action S = J" d 4 xl_(x) + S g f contains the Lagrangian C) written in Eu- 
clidean space and gauge fixing action S g f. Let us consider now a shift: O(x) = 
O b + 0(x) with Ob as a background field and calculate the following expression 
for the partition function in the one-loop approximation: 



Z = 



[DB][Dd][DO+]exp{-S(B,O b ) - 



oO x 



= exp{-F(O b ,g 2 ,^ 2 ,A)}. 
Using the representation 10, we obtain the effective potential: 

Veff = F(4> b) g 2 )H 2 ,A) 



(4) 



(5) 



given by the function F of Eq.10} for the constant background field Ob = c^b = 
const. In this case the one-loop effective potential for monopoles coincides with 
the expression of the effective potential calculated by the authors of Ref.|8| for 
scalar electrodynamics and extended to the massive theory (see review |9[). As it 
was shown in Ref.| 8 |, the effective potential can be improved by consideration of 
the renormalization group equation (RGE). 



2 Renormalization group equations in the Higgs monopole 
model 

The RGE for the effective potential means that the potential cannot depend on a 
change in the arbitrary parameter — renormalization scale M: 

^ = 0. (6) 

The effects of changing it are absorbed into changes in the coupling constants, 
masses and fields, giving so-called running quantities. 

Considering the RG improvement of the effective potential [8,9] and choos- 
ing the evolution variable as 

t = log(cb 2 /M 2 ), (7) 
we have the following RGE for the improved V e ff {$> 2 } with <J) 2 = (J) 2 , |10|: 
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where y is the anomalous dimension and |3 ( ^2 ) , |3 a and |3 g are the RG |3-f unctions 
for mass, scalar and gauge couplings, respectively. RGE JSJ leads to the following 
form of the improved effective potential 1 8 1: 

Veff = -^ un (t)G 2 (t)4> 2 + lA run (t)G 4 (t)4> 4 . (9) 



In our case: 

G(t)=exp[-- 



dt'y(g rutt (t'),A run (t'))]. (10) 



2 _ 

A set of ordinary differential equations (RGE) corresponds to Eq.JSJ: 



dA run 
dt 



(3a (g run (t), Arun(t)) , (11) 



= ^ n (t)P ((tl) ( 8 run(t), A run (t)) , (12) 



run (t), \ruTvW). (13) 

So far as the mathematical structure of HMM is equivalent to HSED, we can use 
all results of the scalar electrodynamics in our calculations, replacing the electric 
charge e and photon field by magnetic charge g and dual gauge field B H . 

The one-loop results for , |3^ 2 ' and y are given in Ref.|8| for scalar field 
with electric charge e, but it is easy to rewrite them for monopoles with charge 
g = gmn: 

T 1 "-^ (14) 



d "ut iI1 ~ = Tfot 2 ^ 39 "" 1 + 10A ™ n ~~ 6Arun 9run). (15) 

^7un ^ p(1] = ^ (4Arun _ 3g 2 un]i (16) 



dt ' ^ 1 16n 2 

^9run ^ o ( 1 ) 9 run /-i y\ 

dt ~ lg ~48tt 2 ' ( ' 

The RG |3-functions for different renormalizable gauge theories with semisim- 
ple group have been calculated in the two-loop approximation and even beyond. 
But in this paper we made use the results of Refs.| 11 1 and 1 12 1 for calculation of 
(3-functions and anomalous dimension in the two-loop approximation, applied 
to the HMM with scalar monopole fields. The higher approximations essentially 
depend on the renormalization scheme. Thus, on the level of two-loop approxi- 
mation we have for all |3-functions: 

(3 = (3' 1) + 8 (2) , (18) 

where 
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and 



The gauge coupling (3 g 2 ' -function is given by Ref. llll : 

n 6 

g(2) _ — y (21) 

Anomalous dimension follows from calculations made in Ref.| 12J: 

y(2) = 3 — n 4 (22) 

7 (16tt 2 ) 2 12 9 ■ ( ' 

In Eqs. il8l — J22I and below, for simplicity, we have used the following notations: 
A = A Tun , g = g run and |X = Hr Un . 



3 The phase diagram in the Higgs monopole model 

Now we want to apply the effective potential calculation as a technique for the 
getting phase diagram information for the condensation of monopoles in HMM. 
If the first local minimum occurs at § — and V e ff (0) = 0, it corresponds to the 
Coulomb-like phase. In the case when the effective potential has the second local 
minimum at cf) — 4> m i n ^ with V^} n (4) 2 rliTl ) < 0, we have the confinement 
phase. The phase transition between the Coulomb-like and confinement phases 
is given by the condition when the first local minimum at <J) = is degenerate 
with the second minimum at cf) = dpo- These degenerate minima are shown in 
Fig.l by the curve 1. They correspond to the different vacua arising in this model. 
And the dashed curve 2 describes the appearance of two minima corresponding 
to the confinement phases. 

The conditions of the existence of degenerate vacua are given by the follow- 
ing equations: 

Veff(0)=Veff(4>g)=0, (23) 



3V e ff, 3V e ff, - ,., ,, 2 , 3V ef f 

-g^-|*=o = -g^-l*=* =0, or V eff (4> ) = — 2-l*=«>o = °- ( 24 ) 
and inequalities 

9 2 Veff , ^ 9^ 

~3$2-'* =0>0, ^ 2_l *=*° ( ' 

The first equation i2"3l applied to Eq.0 gives: 

Hrun =-^run(t )cfoG 2 (to), where t =log(4)§/M 2 ). (26) 
It is easy to find the joint solution of equations 

V eff (ct) 2 )=V^ ff (cl 3 2)=0. (27) 
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Fig.l The effective potential V e g : the curve 1 corresponds 
to the " Coulomb" -" confinement" phase transition; 
curve 2 describes the existence of two minima cor- 
responding to the confinement phases. 



Using RGE {FT}, G3 and Eqs.CSMZZl we obtain: 

Vef f (*o) =^(-Arun(3 ( ^, +A Tun + p A - T A run )G 4 (t )4)o =0, (28) 

or 

(3a +A Tun (l -y- p ((1 . 2) ) =0. (29) 

Substituting in Eq.@9) the functions (3^ 1 , (3 [ 1 \ . and y ( 1 ' given by Eqs.fHl— ATI, 
we obtain in the one-loop approximation the following equation for the phase 
transition border: 

87T 2 

g pt = -2A rurt ( — + A run ) . (30) 

The curve i30l is represented on the phase diagram (A TUTt ; g^ un ) of Fig.2 by the 
curve "1" which describes the border between the Coulomb-like phase with V e f f > 
and the confinement one with V^ n < 0. This border corresponds to the one- 
loop approximation. 

Using Eqs. l|141 - i22^ . we are able to construct the phase transition border in 
the two-loop approximation. Substituting these equations into Eq. l29> . we obtain 
the following phase transition border curve equation in the two-loop approxima- 
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Fig. 2 The one— loop (curve 1) and two— loop (curve 2) ap- 
proximation phase diagram in the dual Abelian Higgs 
model of scalar monopoles. 



tion: 

1 15 97 59 

3u 2 - 1 67t 2 + 6x 2 + ^ (28x 3 + — x 2 y + —xy 2 - -y 3 ) = 0, (31) 

where x = — Apt and y — gp T are the phase transition values of — A run and g 2 un . 
Choosing the physical branch corresponding to g 2 > and g 2 — > 0, when A — > 0, 
we have received curve 2 on the phase diagram (A TU n; g 2 U n.) shown in Fig. 2. This 
curve corresponds to the two-loop approximation and can be compared with 
the curve 1 of Fig.2, which describes the same phase transition border calculated 
in the one-loop approximation. It is easy to see that the accuracy of the 1-loop 
approximation is not excellent and can commit errors of order 30%. 

According to the phase diagram drawn in Fig.2, the confinement phase be- 
gins at g 2 = g^ax and exists under the phase transition border line in the region 
9 2 < 9max' where e 2 is large: e 2 > (27t/g max ) 2 due to the Dirac relation: 

eg=27r, or aa=-^. (32) 



Therefore, we have: 
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in the one-loop approximation: 



9crit 



g 2 max , « 18.61, <W = ^ 



in the two-loop approximation: 



2 

9cTlt 



9max2~15.il, a cr it — "T^T 



1.48, a CTit 



1 .20, a cri t = 



1 



4a, 



crit 



4a c 



0.17 



0.208 
(33) 



Comparing these results, we obtain the accuracy of deviation between them of 
order 20%. 

The last result <!33i coincides with the lattice values obtained for the compact 
QED by Monte Carlo method Q: 



a crtt ~ 0.20 ±0.015, a crit « 1.25 ±0.10. 



(34) 



Writing Eq. <13t with |3 g function given by Eqs.l ll7> . J18L and Pit , we have the 
following RGE for the monopole charge in the two-loop approximation: 



dg 



2 

run 



dt 



9run + 9 run. 



487t 2 (167t 2 ) 2 ' 



or 



dlog a 
dt 



(35) 



The values l3"3l for g 2 Tit = g max i 2 indicate that the contribution of two loops 
described by the second term of Eq. d35l is about 0.3, confirming the validity of 
perturbation theory. 

In general, we are able to estimate the validity of two-loop approximation 
for all (3-functions and y, calculating the corresponding ratios of two-loop con- 
tributions to one-loop contributions at the maxima of curves 1 and 2: 



Acrit — 



9crit 



(2] 



1 



1 ^ 



§2! 



A-"* 1 « -13.16 

9maxl * 18-61 

-0.0080 
-0.0826 

0.1564 



0.3536 



Acrit — 



a 2 • 

"CTlt 



Y 



r 



A max2 w — 7 13 

"run 



9raax2« 15.11 



-0.0065 



-0.0637 



0.0412 



0.2871 



(36) 
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Here we see that all ratios are sufficiently small, i.e. all two-loop contributions 
are small in comparison with one-loop contributions, confirming the validity of 
perturbation theory in the 2-loop approximation. The accuracy of deviation is 
worse (~ 30%) for |3 g -function. But it is necessary to emphasize that calculat- 
ing the border curves 1 and 2 of Fig.2, we have not used RGE J13I for monopole 
charge: |3 g -function is absent in Eq.i2"9l. Therefore, the calculation of gl rit accord- 
ing to Eq. <31> does not depend on the approximation of |3 g function. The above- 
mentioned |3 g -function appears only in the second order derivative of V e ff which 
is related with the monopole mass m (see Refs.[2,3]). 
Eqs.l l33l give the following result: 

«cAt * 5 > (37) 

which is important for the phase transition at the Planck scale predicted by the 
Multiple Point Model (MPM). 

4 Multiple Point Model and Critical Values of the U(l) and 
SU(N) Fine Structure Constants 

Investigating the phase transition in HMM, we had pursued two objects: from 
one side, we had an aim to explain the lattice results, from the other side, we 
were interested in the predictions of MPM. 

4.1 Anti-grand unification theory 

Most efforts to explain the Standard Model (SM) describing well all experimental 
results known today are devoted to Grand Unification Theories (GUTs). The su- 
persymmetric extension of the SM consists of taking the SM and adding the cor- 
responding supersymmetric partners. Unfortunately, at present time experiment 
does not indicate any manifestation of the supersymmetry. In this connection, the 
Anti-Grand Unification Theory (AGUT) was developed in Refs. [13-17, 4] as a re- 
alistic alternative to SUSY GUTs. According to this theory, supersymmetry does 
not come into the existence up to the Planck energy scale: Mp\ = 1 .22 ■ 1 1 9 GeV. 
The Standard Model (SM) is based on the group SMG: 

SMG = SU(3) C x SU(2)l x U(1)y- (38) 

AGUT suggests that at the energy scale |Xg ~ M-pv = M-pi there exists the more fun- 
damental group G containing lM gen copies of the Standard Model Group SMG: 

G = SMG, x SMG 2 x ... x SMG Ngen = (SMG) N «« n , (39) 

where N gen designates the number of quark and lepton generations (families). 
If N g en =3 (as AGUT predicts), then the fundamental gauge group G is: 

G = (SMG) 3 = SMGi st gen. X SMG 2 ndgen. X SMG 3 rdgen., (40) 
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or the generalized ones: 

G f = (SMG) 3 x U(1) f , or G ext =(SMGxU(l) B _ L ) 3 , (41) 

which were suggested by the fitting of fermion masses of the SM (see Refs.|16|), 
or by the see-saw mechanism with right-handed neutrinos 1 17 \. 

4.2 Multiple Point Principle 

AGUT approach is used in conjuction with the Multiple Point Principle proposed 
in Ref .|4J. According to this principle Nature seeks a special point — the Multiple 
Critical Point (MCP) — which is a point on the phase diagram of the fundamental 
regulirized gauge theory G (or Gf, or G ex t), where the vacua of all fields existing 
in Nature are degenerate having the same vacuum energy density Such a phase 
diagram has axes given by all coupling constants considered in theory. Then all 
(or just many) numbers of phases meet at the MCP. MPM assumes the existence 
of MCP at the Planck scale, insofar as gravity may be "critical" at the Planck scale. 

The philosophy of MPM leads to the necessity to investigate the phase tran- 
sition in different gauge theories. A lattice model of gauge theories is the most 
convenient formalism for the realization of the MPM ideas. As it was mentioned 
above, in the simplest case we can imagine our space-time as a regular hyper- 
cubic (3+l)-lattice with the parameter a equal to the fundamental (Planck) scale: 
a = A P = 1 /Mpi. 



4.3 AGUT-MPM prediction of the Planck scale values of the U(l), SU(2) and 
SU(3) fine structure constants 

The usual definition of the SM coupling constants: 

5 OC CC Q 2 

3 cos 2 8^ sin 2 4tt 

where a and a s are the electromagnetic and strong fine structure constants, re- 
spectively, is given in the Modified minimal subtraction scheme (MS). Here 9^rs~ 
is the Weinberg weak angle in MS scheme. Using RGE with experimentally es- 
tablished parameters, it is possible to extrapolate the experimental values of three 
inverse running constants oc^ 1 ([i] (here [i is an energy scale and i=l,2,3 corre- 
spond to U(l), SU(2) and SU(3) groups of the SM) from the Electroweak scale to 
the Planck scale. The precision of the LEP data allows to make this extrapolation 
with small errors (see |18|). Assuming that these RGEs for a^ 1 contain only 
the contributions of the SM particles up to (a. w u. Pl and doing the extrapolation 
with one Higgs doublet under the assumption of a "desert", the following results 
for the inverses ocy \ 3 (here ccy = §oci ) were obtained in Ref. ] 4 1 (compare with 
EH): 

ocyVpi) «55.5; 1 (hpi) « 49.5; 1 (ji Pl ) « 54.0. (43) 
The extrapolation of oiy 2 3 (u.) up to the point |x = [ip\ is shown in Fig. 3. 
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Fig. 3 The evolution of three inverse running constants a y 1 (fi), 
where i=l,2,3 correspond to U(l), SU(2) and SU(3) groups 
of the SM. The extrapolation of their experimental values 
from the Electroweak scale to the Planck scale was obtained 
by using the renormlization group equations with one Higgs 
doublet under the assumption of a "desert". The precision 
of the LEP data allows to make this extrapolation with 
small errors. AGUT works in the region (j, G < /x < jip/. 



According to AGUT, at some point \i — \iq < \in (but near \in) the funda- 
mental group G (or Gf, or G ex ^) undergoes spontaneous breakdown to its diago- 
nal subgroup: 

G — > Gdiag.subgr. = {g, g, 9 1 1 9 £ SMG}, (44) 

which is identified with the usual (low-energy) group SMG. 

The AGUT prediction of the values of cti{\±) at |j. = upi is based on the MPM 
assumptions, and gives these values in terms of the corresponding critical cou- 
plings a^crit [13-15,4]: 

OilH-Pl) = -TT = for 1 = 2 > 3 > ( 45 ) 

N gen J 

and 

ai(^Pi) = y- — — TT = — 7— for U(1). (46) 

lN gen (N gen + 1) 6 

It was assumed in Ref.|4| that the MCP values ai iCr it in Eqs.l45l and j46l coincide 
with the triple point values of the effective fine structure constants given by the 
lattice SU(3)-, SU(2)- and U(l)-gauge theories. 
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If the point \i = \ig is very close to the Planck scale [i — [ipi, then according 
to Eqs.133 and l4"6l . we have: 

a lst gen. ~ a 2nd gen. ~ a 3rd gen. ~ /■ ~ 7 > V*' / 



what is almost equal to the value: 



(48) 



obtained theoretically by Parisi improvement method for the Coulomb-like phase 
[4,6]. The critical value (48) is close to the lattice and HMM ones: see Eq. <37> . This 
means that in the U(l) sector of AGUT we have a near the critical point, and we 
can expect the existence of MCP at the Planck scale. 
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Abstract. We review our recent development of family replicated gauge group model, 
which generates the Large Mixing Angle MSW solution. The model is based on each family 
of quarks and leptons having its own set of gauge fields, each containing a replica of the 
Standard Model gauge fields plus a (B — L)-coupled gauge field. A fit of all the seventeen 
quark-lepton mass and mixing angle observables, using just six new Higgs field vacuum 
expectation values, agrees with the experimental data order of magnitudewise. However, 
this model can not predict the baryogenesis in right order, therefore, we discuss further 
modification of our model and present a preliminary result of baryon number to entropy 
ratio. 

1 Introduction 

We have previously attempted to fit all the fermion masses and their mixing an- 
gles 1 1 2 1 including baryogenesis 1 3 1 in a model without supersymmetry or grand 
unification. This model has the maximum number of gauge fields consistent with 
maintaining the irreduciblity of the usual Standard Model fermion representa- 
tions, added three right-handed neutrinos. The predictions of this previous model 
are in order of magnitude agreement with all existing experimental data, how- 
ever, only provided we use the Small Mixing Angle MSW |4| (SMA-MSW) solu- 
tion. But, for the reasons given below, the SMA-MSW solution is now disfavoured 
by experiments. So here we review a modified version of the previous model, 
which manages to accommodate the Large Mixing Angle MSW (LMA-MSW) so- 
lution for solar neutrino oscillations using 6 additional Higgs fields (relative to 
the Standard Model) vacuum expectation values (VEVs) as adjustable parame- 
ters. 

A neutrino oscillation solution to the solar neutrino problem and a favouring 
of the LMA-MSW solution is supported by SNO results 1 5 1 : The measurement of 
the 8 B and hep solar neutrino fluxes shows no significant energy dependence 
of the electron neutrino survival probability in the Super-Kamiokande and SNO 
energy ranges. 

Moreover, the important result which also supports LMA-MSW solution on 
the solar neutrino problem, reported by the Super-Kamiokande collaboration 1 6 1, 

* E-mail: hbech@mail.desy.de 
** E-mail: yasutaka@mail.desy.de 
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that the day-night asymmetry data disfavour the SMA-MSW solution at the 95% 
C.L.. 

In fact, global analyses I7I8I9I10I of all solar neutrino data have confirmed 
that the LMA-MSW solution gives the best fit to the data and that the SMA-MSW 
solution is very strongly disfavoured and only acceptable at the 3o~ level. Typical 
best fit values of the mass squared difference and mixing angle parameters in the 
two flavour LMA-MSW solution are Am| w 4.5 x 10" 5 eV 2 and tan 2 9© w 0.35. 

This paper is organised as follows: In the next section, we present our gauge 
group - the family replicated gauge group - and the quantum numbers of fermion 
and Higgs fields. Then, in section 3 we discuss our philosophy of all gauge- and 
Yukawa couplings at Planck scale being of order unity. In section 4 we address 
how the family replicated gauge group breaks down to Standard Model gauge 
group, and we add a small review of see-saw mechanism. The mass matrices of 
all sectors are presented in section 5, the renormalisation group equations - renor- 
malisable and also 5 dimensional non-renormalisable ones - are shown in section 
6. The calculation is described in section 7 and the results are presented in sec- 
tion 8. We discuss further modification of our model and present a preliminary 
results of baryon number to entropy ratio in section 9. Finally, section 1 contains 
our conclusion. 

2 Quantum numbers of model 

Our model has, as its back-bone, the property that there are generations (or fami- 
lies) not only for fermions but also for the gauge bosons, i.e., we have a generation 
(family) replicated gauge group namely 

^^(SMGixUO),^) , (1) 

where SMG denotes the Standard Model gauge group = SU(3) x SU(2) x U(l ), 
x denotes the Cartesian product and i runs through the generations, i = 1,2,3. 

Note that this family replicated gauge group, eq. Q, is the maximal gauge 
group under the following assumptions: 

• It should only contain transformations which change the known 45 (= 3 gen- 
erations of 15 Weyl particles each) Weyl fermions of the Standard Model 
and the additional three heavy see-saw (right-handed) neutrinos. That is our 
gauge group is assumed to be a subgroup of U(48). 

• We avoid any new gauge transformation that would transform a Weyl state 
from one irreducible representation of the Standard Model group into another 
irreducible representation: there is no gauge coupling unification. 

• The gauge group does not contain any anomalies in the gauge symmetry - 
neither gauge nor mixed anomalies even without using the Green-Schwarz 
anomaly cancelation mechanism. 

• It should be as big as possible under the foregoing assumptions. 

The quantum numbers of the particles/ fields in our model are found in table 
1 and use of the following procedure: In table 1 one finds the charges under the 
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Table 1. All U(1) quantum charges in the family replicated model. The symbols for the 
fermions shall be considered to mean "proto"-particles. Non-abelian representations are 
given by a rule from the abelian ones (see Eq. 
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six U(1) groups in the gauge group Then for each particle one should take 
the representation under the SU(2)i and SU(3)i groups (i = 1,2,3) with lowest 
dimension matching to \ji/2 according to the requirement 

T + + ¥ = (modl) ' (2) 

where t, and dt are the triality and duality for the i'th proto-generation gauge 
groups SU(3)t and SU(2)i respectively. 

3 The philosophy of all couplings being order unity 

Any realistic model and at least certainly our model tends to get far more fun- 
damental couplings than we have parameters in the Standard Model and thus 
pieces of data to fit. This is especially so for our model based on many U(1) 
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charges |11| because we take it to have practically any not mass protected par- 
ticles one may propose at the fundamental mass scale, taken to be the Planck 
mass. Especially we assume the existence of Dirac fermions with order of funda- 
mental scale masses needed to allow the quark and lepton Weyl particles to take 
up successively gauge charges from the Higgs fields VEVs. So unless we make as- 
sumptions about the many coupling constants and fundamental masses we have 
no chance to predict anything. Almost the only chance of making an assumption 
about all these couplings, which is not very model dependent, is to assume that 
they are all of order unity in the fundamental unit. This is the same type of as- 
sumption that is really behind use of dimensional arguments to estimate sizes of 
quantities. A procedure very often used successfully. If we really assumed every 
coupling and mass of order unity we would get the effective Yukawa couplings 
of the quarks and leptons to the Weinberg-Salam Higgs field to be also of or- 
der unity what is phenomenologically not true. To avoid this prediction we then 
blame the smallness of all but the top-Yukawa coupling on smallness in funda- 
mental Higgs VEVs. That is to say we assume that the VEVs of the Higgs fields in 
Table 1, p, cu, T, W, x, 4>b-l an d 4>ws ar e (possibly) very small compared to the 
fundamental /Planck unit, and these are the quantities we have to fit. 

Technically we implement these unknown - but of order unity according to 
our assumption - couplings and masses by a Monte Carlo technique: we put them 
equal to random numbers with a distribution dominated by numbers of order 
unity and then perform the calculation of the observable quantities such as quark 
or lepton masses and mixing angles again and again. At the end we average the 
logarithmic of these quantities and exponentiate them. In this way we expect 
to get the typical order of magnitude predicted for the observable quantities. In 
praxis we do not have to put random numbers in for all the many couplings in 
the fundamental model, but can instead just provide each mass matrix element 
with a single random number factors. 

After all a product of several of order unity factors is just an order unity 
factor again. To resume our model philosophy: Only Higgs field VEVs are not of 
order unity. We must be satisfied with order of magnitude results. 

4 Breaking of the family replicated gauge group to the Standard 
Model 

The family replicated gauge group broken down to its diagonal subgroup at 
scales about one or half order of magnitude under the Planck scale by Higgs 
fields - W, T, cu, p and x (in Table^: 

t _* 2 , 3 (SMGi x UCIJb^.O -» SMG x U(1)b-l ■ (3) 

This diagonal subgroup is further broken down by yet two more Higgs fields — 
the Weinberg-Salam Higgs field 4> W s and another Higgs field 4>b-l — to SU(3) x 
U(l) em . 



146 H. B. Nielsen and Y. Takanishi 



4.1 See-saw mechanism 

See-saw mechanics is build into our model to fit the scale of the neutrino oscilla- 
tions, i.e., we use the right-handed neutrinos with heavy Majorana masses (10 1 1 
GeV). 

In order to mass-protect the right-handed neutrino from getting Planck scale 
masses, we have to introduce cf) B _ L which breaks the B — L quantum charge spon- 
taneously, and using this new Higgs filed we are able to deal the neutrino oscil- 
lations, i.e., to fit the scale of the see-saw particle masses. However, due to mass- 
protection by the Standard Model gauge symmetry, the left-handed Majorana 
mass terms should be negligible in our model. Then, naturally, the light neutrino 
mass matrix - effective left-left transition Majorana mass matrix - can be obtained 
via the see-saw mechanism |12|: 

Meff^M^MR 1 (M°) T . (4) 



5 Mass matrices 

Using the U(l ) fermion quantum charges and Higgs field (presented in Table 
we can calculate the degrees of suppressions of the left-right transition - Dirac 
mass - matrices and also Majorana mass matrix (right-right transition). 

Note that the random complex order of unity numbers which are supposed 
to multiply all the mass matrix elements are not represented in following matri- 
ces: the up-type quarks: 

(cut) 3 WtT 2 tuptW^T 2 tupt(Wt) 2 T^ 

)4 pW tT 2 WtT 2 (Wt) 2 T j (5) 



cu 3 W(Tt) 2 tVW(Tt) 2 tup t T 3N 

pW(Tt) 2 W(Tt) 2 T 3 I (6) 




w 2 pW 2 {V) 4 W 2 (T 




the charged leptons: 



, , / u) 3 W(Tt) 2 (cu t ) 3 p 3 W(Tt) 2 (a)t) 3 p 3 WT 4 x\ 
M E ~iM| u) 6(ptj3 W(T t)2 w(Tt) 2 wr*x I (7) 

t)3 (w t)2 T 4 ( W t) 2 T 4 Wi- 




the Dirac neutrinos: 




tut) 3 p 3 WtT 2 (cu t ) 3 p 3 WtT 2 x\ 

WtT 2 W^T 2 x (8) 

WtV W+Tt / 



and the Majorana (right-handed) neutrinos: 




(9) 
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6 Renormalisation group equations from Planck scale to week 
scale via see-saw scale 

It should be kept in mind that the effective Yukawa couplings for the Weinberg- 
Salam Higgs field, which are given by the Higgs field factors in the above mass 
matrices multiplied by order unity factors, are the running Yukawa couplings 
at a scale near the Planck scale. In this way, we had to use the renormalisation 
group (one-loop) P-functions to run these couplings down to the experimentally 
observable scale which we took for the charged fermion masses to be compared to 
"measurements" at the scale of 1 GeV, except for the top quark mass prediction. 
We define the top quark pole mass: 



where we put M = 1 80 GeV for simplicity. 

We use the one-loop (3 functions for the gauge couplings and the charged 
fermion Yukawa matrices |13| as follows: 



where t = In \i. 

By calculation we use the following initial values of gauge coupling con- 
stants: 




(10) 




(11) 



Tr(3Yt Y u +3Yt Y D +YtY E ) , 



U(l) 
SU(2) 
SU(3) 



91 (M z ) = 0.462 , 9l (M Planck ) = 0.614 
gi(Mz) = 0.651 , g 2 (M Planck ) = 0.504 
g 3 (M z ) = 1.22 , g 3 (M Planck ) =0.491 



(12) 
(13) 
(14) 
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6.1 The renormalisation group equations for the effective neutrino mass 
matrix 

The effective light neutrino masses are given by an irrelevant, nonrenormal- 
isable (5 dimensional term) - effective mass matrix M e ff - for which the running 
formula is [14 J: 

167t2 ^T = (-392 + 2A + 2Y s )M eff -^ (M eff (Y E Yt) T + (Y E Yt) M eff ) , (15) 

where A is the Weinberg-Salam Higgs self -coupling constant and the mass of the 
Standard Model Higgs boson is given by M^ = A (4>ws) 2 - We just for simplicity 
take Mh = 115 GeV thereby we ignore the running of the Higgs self-coupling 
and fixed as A = 0.2185. 

Note that the renormalisation group equations are used to evolve the effec- 
tive neutrino mass matrix from the see-saw sale, set by (4)b-l) in our model, to 
1 GeV. 



7 Method of numerical computation 

In the philosophy of order unity numbers spelled out in section [3] we evaluate 
the product of mass-protecting Higgs VEVs required for each mass matrix ele- 
ment and provide it with a random complex number, Atj, of order one as a factor 
taken to have Gaussian distribution with mean value zero. But we hope the exact 
form of distribution does not matter much provided we have (In |A t j |) = 0. In 
this way, we simulate a long chain of fundamental Yukawa couplings and prop- 
agators making the transition corresponding to an effective Yukawa coupling in 
the Standard Model and the parameters in neutrino sector. In the numerical com- 
putation we then calculate the masses and mixing angles time after time, using 
different sets of random numbers and, in the end, we take the logarithmic average 
of the calculated quantities according to the following formula: 



™H ex P 1^ • (16) 




Here (m) is what we take to be the prediction for one of the masses or mixing 
angles, mi is the result of the calculation done with one set of random number 
combinations and N is the total number of random number combinations used. 

Since we only expect to make order of magnitude fits, we should of course 
not use ordinary x 2 defined form the experimental uncertainties by rather the 
X 2 that would correspond to a relative uncertainly - an uncertain factor of order 
unity. Since the normalisation of such a x 2 is not so easy to choose exactly we de- 
fine instead a quantity which we call the goodness of fit (g.o.f.). Since our model 
can only make predictions order of magnitudewise, this quantity g.o.f. should 
only depend on the ratios of the fitted masses and mixing angles to the experi- 
mentally determined masses and mixing angles: 



;.o.f. ^ Y_ 



In 



(m) 



(17) 
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where m eX p are the corresponding experimental values presented in Table|2 

We should emphasise that we do not adjust the order of one numbers by se- 
lection, i.e., the complex random numbers are needed for only calculational pur- 
poses. That means that we have only six adjustable parameters - VEVs of Higgs 
fields - and, on the other hand, that the averages of the predicted quantities, (m), 
are just results of integration over the "dummy" variables - random numbers - 
therefore, the random numbers are not at all parameters! 

Strictly speaking, however, one could consider the choice of the distribution 
of the random order unity numbers as parameters. But we hope that provided we 
impose on the distribution the conditions that the average be zero and the aver- 
age of the logarithm of the numerical value be zero, too, any reasonably smooth 
distribution would give similar results for the (m) values at the end. In our early 
work [2J we did see that a couple of different proposals did not make too much 
difference. 

8 Results 

We averaged over N = 10, 000 complex order unity random number combina- 
tions. These complex numbers are chosen to be a number picked from a Gaussian 
distribution, with mean value zero and standard deviation one, multiplied by a 
random phase factor. We put them as factors into the mass matrices (eqs. 15191 . 
Then we computed averages according to eq. J16I and used eq. \Y7\ as a x 2 to fit 
the 6 free parameters and found: 

(4) ws ) = 246GeV, (<|> B _ L ) =1.64x 10 11 GeV, (u>) = 0.233, 

(p) = 0.246, (W) =0.134, (T) = 0.0758, ( x ) =0.0737, (18) 

where, except for the Weinberg-Salam Higgs field and (<t> B _ L ), the VEVs are ex- 
pressed in Planck units. Hereby we have considered that the Weinberg-Salam 
Higgs field VEV is already fixed by the Fermi constant. The results of the best fit, 
with the VEVs in eq. | |18> , are shown in Tableland the fit has g.o.f. = 3.63. 

We have 11 = 17 — 6 degrees of freedom - predictions - leaving each of 
them with a logarithmic error of ^3.63/11 ~ 0.57, which is very close to the 
theoretically expected value 64% |15|. This means that we can fit all quantities 
within a factor exp f^/3.63/11^ ~ 1 .78 of the experimental values. 

From the table |2 the experimental mass values are a factor two higher than 
predicted for down, charm and for the Cabibbo angle V us while they are smaller 
by a factor for strange and electron. Thinking only on the angles and masses (not 
squared) the agreement is in other cases better than a factor two. 

Experimental results for the values of neutrino mixing angles are often pre- 
sented in terms of the function sin 2 29 rather than tan 2 (which, contrary to 
sin 2 20, does not have a maximum at = 7t/4 and thus still varies in this re- 
gion). Transforming from tan 2 variables to sin 2 20 variables, our predictions for 
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Table 2. Best fit to conventional experimental data. All masses are running masses at 1 GeV 
except the top quark mass which is the pole mass. Note that we use the square roots of the 
neutrino data in this Table, as the fitted neutrino mass and mixing parameters (m) , in our 
goodness of fit (g.o.f.) definition, eq. 1171 . 
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the neutrino mixing angles become: 

sin 2 26 Q = 0.61 , (19) 
sin 2 26 atm =0.96, (20) 
sin 2 20,3 =0.17. (21) 

We also give here our predicted hierarchical neutrino mass spectrum: 

(22) 
(23) 
(24) 

Our agreement with experiment is excellent: all of our order of magnitude 
neutrino predictions lie inside the 99% C.L. border determined from phenomeno- 
logical fits to the neutrino data, even including the CHOOZ upper bound. Our 
prediction of the solar mass squared difference is about a factor of 2 larger than 
the global data fit even though the prediction is inside of the LMA-MSW region, 
giving a contribution to our goodness of fit of g.o.f. w 0.14. Our CHOOZ angle 
also turns out to be about a factor of 2 larger than the experimental limit at 90% 
C.L., delivering another contribution of g.o.f. w 0.14. In summary our predictions 
for the neutrino sector agree extremely well with the data, giving a contribution 
of only 0.34 to g.o.f. while the charged fermion sector contributes 3.29 to g.o.f.. 
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8.1 CP violation 

Since we have taken our random couplings to be - whenever allowed - complex 
we have order of unity or essentially maximal CP-violation so a unitary triangle 
with angles of order one is a success of our model. After our fitting of masses and 
of mixings we can simply predict order of magnitudewise of CP-violation in e.g. 
K° — K° decay or in CKM and MNS mixing matrices in general. 

The Jarlskog area J CP provides a measure of the amount of CP violation in 
the quark sector 1 16 [ and, in the approximation of setting cosines of mixing angles 
to unity, is just twice the area of the unitarity triangle: 

Jcp = V us V cb V ub sin 6 , (25) 

where 6 is the CP violation phase in the CKM matrix. In our model the quark 
mass matrix elements have random phases, so we expect 6 (and also the three 
angles a, (3 and y of the unitarity triangle) to be of order unity and, taking an 
average value of | sin 6 ~ 1 /2, the area of the unitarity triangle becomes 

Jcp « ^V us V cb V ub . (26) 

Using the best fit values for the CKM elements from Table El we predict J CP w 
3.1 x 10~ 6 to be compared with the experimental value (2 — 3.5) x 10~ 5 . Since 
our result for the Jarlskog area is the product of four quantities, we do not expect 
the usual ±64% logarithmic uncertainty but rather ±V4~ ■ 64% = 128% logarith- 
mic uncertainty. This means our result deviates from the experimental value by 
^ 3 ixio-^ /1 -28 = 1.7 "standard deviations". 

The Jarlskog area has been calculated from the best fit parameters in Table 
2, it is also possible to calculate them directly while making the fit. So we have 
calculated J CP for 1M = 10,000 complex order unity random number combina- 
tions. Then we took the logarithmic average of these 10, 000 samples of J CP and 
obtained the following result: 

Jcp =3.1 x 10- 6 , (27) 
in good agreement with the values given above. 



8.2 Neutrinoless double beta decay 



Another prediction, which can also be made from this model, is the electron 
"effective Majorana mass" - the parameter in neutrinoless beta decay - defined 
by: 



|(m) 



TTV; 



(28) 



where mi are the masses of the neutrinos "Vt and U e i are the MNS mixing ma- 
trix elements for the electron flavour to the mass eigenstates i. We can substitute 
values for the neutrino masses mi from eqs. ( 1221241 and for the fitted neutrino 
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mixing angles from Table |3 into the left hand side of eq. t28> . As already men- 
tioned, the CP violating phases in the MNS mixing matrix are essentially random 
in our model. So we combine the three terms in eq. i28\ by taking the square root 
of the sum of the modulus squared of each term, which gives our prediction: 



In the same way as being calculated the Jarlskog area we can compute using 
N = 10, 000 complex order unity random number combinations to get the |(m)|. 
Then we took the logarithmic average of these 1 0, 000 samples of | (m) | as usual: 



This result does not agree with the central value of recent result - "evidence" - 
from the Heidelberg-Moscow collaboration \17\. 

9 Baryogenesis via Lepton Number Violation 

Having now a well fitted model giving orders of magnitude for all the Yukawa 
couplings and having the see-saw mechanism, it is obvious that we ought to 
calculate the amount of baryons Yb relative to entropy being produced via the 
Fukugita-Yanagida mechanism |18|. According to this mechanism the decay of 
the right-handed neutrinos by CP-violating couplings lead to an excess of the 
B — L charge (meaning baryon number minus lepton number), the relative ex- 
cess in the decay from Majorana neutrino generation number i being called et. 
This excess is then immediately - and continuously back and forth - being con- 
verted partially to a baryon number excess, although it starts out as being a 
lepton number L asymmetry, since the right-handed neutrinos decay to leptons 
and Weinberg-Salam Higgs particles. It is a complicated discussion to estimate to 
what extend the B — L asymmetry is washed out later in the cosmological devel- 
opment, but our estimates goes that there is not enough baryon number excess 
left to fit the Big Bang development at the stage of formation of the light elements 
primordially (nuclearsyn thesis). 

Recently we have, however, developed a modified version |19| of our model 
- only deviating in the right-handed sector - characterized by changing the quan- 
tum numbers assumed for the see-saw scale producing Higgs field 4>b-l in such 
a way that the biggest matrix elements in the right-handed mass matrix (eq. |5} 
becomes the pair of - because of the symmetry - identical off diagonal elements 
(row, column)=(2,3) and (3,2). Thereby we obtain two almost degenerate right- 
handed neutrinos and that helps for making the B — L asymmetry in the decay 
bigger. In this modified model that turns out to fit the rest of our predictions ap- 
proximately equally well or even better we then get a very satisfactory baryon 
number relative to entropy prediction 



|(m)|«3.1 x 10~ 3 eV. 



(29) 



|(m)| =4.4 x 10~ 3 eV. 



(30) 



Y B 



2.5 x 10 



1 1 



(31) 



In the same time as making this modification of the 4>b-l quantum numbers 
we also made some improvements in the calculation by taking into account the 
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running of the Dirac neutrino Yukawa couplings from the Planck scale to the 
corresponding right-handed neutrino scales. Also, we calculated more accurate 
dilution factors than previous our work |3 [. However, foregoing work was based 
on the mass matrices which predicted the SMA-MSW, so we must investigate 
the baryogenesis using the present mass matrices, of course, with the modified 
right-handed Majorana mass matrix. 

10 Conclusion 

We have reviewed our model which is able to predict the experimentally fa- 
vored LMA-MSW solution rather than the SMA-MSW solution for solar neutrino 
oscillations after careful choice of the U(l) charges for the Higgs fields causing 
transitions between 1st and 2nd generations. However, the fits of charged lepton 
quantities become worse compare to our "old" model that can predict SMA-MSW 
solar neutrino solution. On the other hand, we now can fit the neutrino quanti- 
ties very well: the price paid for the greatly improved neutrino mass matrix fit - 
the neutrino parameters now contribute only very little to the g.o.f. - is a slight 
deterioration in the fit to the charged fermion mass matrices. In particular the 
predicted values of the quark masses ma and m c and the Cabibbo angle V us are 
reduced compared to our previous fits. However the overall fit agrees with the 
seventeen measured quark-lepton mass and mixing angle parameters in Table 2 
within the theoretically expected uncertainty |15[ of about 64%; it is a perfect fit 
order of magnitudewise. It should be remarked that our model provides an order 
of magnitude fit /understanding of all the effective Yukawa couplings of the Stan- 
dard Model and the neutrino oscillation parameters in terms of only 6 parameters 
- the Higgs field vacuum expectation values! 
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and The Niels Bohr Institute, Blegdamsvej 17, Copenhagen 0, Denmark 

Abstract. In our model with a Standard Model gauge group extended with a baryon 
number minus lepton number charge for each family of quarks and leptons, we calculate the 
baryon number relative to entropy produced in early Big Bang by the Fukugita-Yanagida 
mechanism. With the parameters, i.e., the Higgs VEVs already fitted in a very success- 
ful way to quark and lepton masses and mixing angles we obtain the order of magnitude 
pure prediction Y B =2.59 _2i5 x 10 which according to a theoretical estimate should 
mean in this case an uncertainty of the order of a factor 7 up or down (to be compared to 
Y B = (1.7-8.1) x 10~") using a relatively crude approximation for the dilution factor, 
while using another estimate based on Buchmiiller and Pliimacher a factor 500 less, but 
this should rather be considered a lower limit. With a realistic uncertainty due to wash-out 
of a factor 1 00 up or down we even with the low estimate only deviate by 1 .5a. 

1 Introduction 

Using the model for mass matrices presented by us in an other contribution 1 1 1 
at this conference we want to compute the amount of baryons produced in the 
early universe. This model works by having the mass matrix elements being sup- 
pressed by approximately conserved quantum numbers from a gauge group re- 
peated for each family of quarks and leptons and also having a (B — L) charge for 
each family. 

The baryon number density relative to entropy density, Yb, is one of the 
rather few quantities that can give us information about the laws of nature be- 
yond the Standard Model and luckily we have from the understanding of the 
production of light isotopes at the minute scale in Big Bang fits to this quan- 
tity 1 2 1 . The "experimental" data of the ratio of baryon number density to the 
entropy density is 



We already had a good fit of all the masses and mixings 1 2 1 [ for both quarks and 
leptons measured so far and agreeing with all the bounds such as neutrinoless 
beta decay and proton decay not being seen and matching on the borderline but 
consistent with the accuracy of our model and of the experiment of CHOOZ the 
electron to heaviest left-handed neutrino mixing, and that in a version of our 

* E-mail: hbech@mail.desy.de 
** E-mail: yasutaka@mail.desy.de 



Y B 



(1.7-8.1) x 10 



1 1 



(1) 



exp 
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model in which the dominant matrix element in the right-handed neutrino mass 
matrix is the diagonal one for the "third" (i.e. with same (B — L)i as the third 
family) family -vr3 right-handed neutrino. This version of our model which fits 
otherwise very well does not give sufficient (B — L) excess, that survives, but the 
by now the best model in our series should have the right-handed mass matrix 
dominated by the off-diagonal elements (2,3) and (3,2), so that there appears 
two almost mass degenerate see-saw neutrinos, in addition to the third one (first 
family) which is much lighter. 

2 Mass matrices and results for masses and mixing angles 

Our model produces mass matrix elements - or effective Yukawa couplings - 
which are suppressed from being of the order of the top-mass because they are 
forbidden by the conservation of the gauge charges of our model and can only be- 
come different from zero using the 6 Higgs fields 1 1 3 [ which we have in addition 
to the field replacing the Weinberg-Salam one. In the neutrino sector according to 
the see-saw mechanism 1 6 [ we have to calculate Dirac- and Majorana-mass ma- 
trices, MeffW M° Mj, 1 (M° ) T , to obtain the effective mass matrix M eff for the left 
handed neutrinos we in practice can "see". Here we present all mass matrices as 
they follow from our choice of quantum numbers for the 7 Higgs fields in our 
model and for the quarks and leptons (as they can be found in the other contri- 
bution). Only the quantum numbers for the field called 4>b-l is - in order to get 
degenerate see-saw neutrinos - changed into having the B — L quantum numbers 
of family 2 and 3 equal to 1, i.e., (B — L)2 = (B — LJ3 = 1, while the other family 
quantum numbers are just zero: 
the up-type quarks: 




o>t) 4 pWtT 2 Wn 2 (Wt) 2 T 
(cut) 4 p 1 wtTt 



(2) 



the down- type quarks: 




w 3 W{V) 2 cuptW(Tt) 2 cuptT 3 
v 2 pW{V} 2 W(V) 2 T 3 
u 2 pW 2 (Tt) 4 W 2 (Tt) 4 WT 



(3) 



the charged leptons: 




(cut) 3 p 3 W(Tt) 2 (cutjVwrY 
W(Tt) 2 WT 4 X 
(W 1 ') 2 !" 4 WT 



(4) 



the Dirac neutrinos: 




(cu t ) 3 WtT 2 (tut) 3 p 3 WtT 2 (cu t ) 3 p 3 WtT 2 x x 
(pt)3 W t T 2 w t T 2 WtT 2 x 
pt)3 w t T t x t wtTV W+Tt , 



(5) 



Family Replicated Calculation of Baryogenesis 



157 



and the Majorana (right-handed) neutrinos: 



/ (pW (pW/2 (p^ 3 /^ 
Mr~(c|> b _ l > (pW/2 X 1 1 

V (P f ) 3 /2 1 X , 



(6) 



We shall remember that it is here understood that all the matrix elements are to be 
provided with order of unity factors which we do not know and in practice have 
treated by inserting random order of unity factors over which we then average at 
the end (in a logarithmic way). 

3 Renormalisation group equations 

The model for the Yukawa couplings we use gives, in principle, these couplings 
at the fundamental scale, taken to be the Planck scale, at first, and we then use 
the renormalisation group to run them down to the scales where they are to be 
confronted with experiment. From the Planck scale down to the see-saw scale or 
rather from where our gauge group is broken down to SMG x U(1)b-l we use 
the one-loop renormalisation group running of the Yukawa coupling constant 
matrices and the gauge couplings 1 3 1 in GUT notation including the running of 
Dirac neutrino Yukawa coupling: 



where t = In u. and |J. is the renormalisation point. 

In order to run the renormalisation group equations down to 1 GeV, we use 
the following initial values: 



We varied the 6 free parameters and found the best fit, corresponding to the low- 




U ) 



U(1) 
SU(2) 
SU(3) 



gi(M z )= 0.462 , gi (M Planck ) = 0.614 , 
g 2 (M z ) = 0.651 , g 2 (M Planck ) =0.504, 
g 3 (M z ) = 1.22 , g 3 (M Planck ) = 0.491. 



(7) 
(8) 
(9) 




n 2 



3.38, with the following 



(4> ws ) =246 GeV , (4> B _ L > = 1 .23 x 10 10 GeV , (o>) = 0.245 , 
(p) = 0.256 , (W) = 0.143 , (T) = 0.0742 , (x) = 0.0408 , (10) 
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where, except for the Weinberg-Salam Higgs field and (4> B _ t ), the VEVs are ex- 
pressed in Planck units. Hereby we have considered that the Weinberg-Salam 
Higgs field VEV is already fixed by the Fermi constant. The results of the best fit, 
with the VEVs in eq. < fT0i . are shown in Tabled 

4 Quantities to use for baryogenesis calculation 

Since the baryogenesis in the Fukugita-Yanagida scheme [7\ arises from a nega- 
tive excess of lepton number being converted by Sphalerons to a positive baryon 
number excess partly and this negative excess comes from the CP violating decay 
of the see-saw neutrinos we shall introduce the parameters et giving the mea- 
sure of the relative asymmetry under C or CP in the decay of neutrino number i: 
Defining the measure £i for the CP violation 

£l L«, p r(N I , l ->£«(|)P s ) + 5: a ,pr(N llt ->l*(|>gt) ' 

where V are N R i decay rates (in the N „ t rest frame), summed over the neutral and 
charged leptons (and Weinberg-Salam Higgs fields) which appear as final states 
in the N R i decays one sees that the excess of leptons over anti-leptons produced 





Fitted 


Experimental 


m u 


5.2 MeV 


4 MeV 


m d 


5.0 MeV 


9 MeV 


m e 


1.1 MeV 


0.5 MeV 


m c 


0.70 GeV 


1 .4 GeV 


m s 


340 MeV 


200 MeV 


nV 


81 MeV 


105 MeV 


M t 


208 GeV 


180 GeV 


TRb 


7.4 GeV 


6.3 GeV 


TUt 


1.11 GeV 


1.78 GeV 


Vus 


0.10 


0.22 


V cb 


0.024 


0.041 


Vub 


0.0025 


0.0035 


Am 


9.0 x lO - * eV 2 


4.5 x 10~ s eV 2 


AmL 


1.8 x 10~ 3 eV 2 


3.0 x 10~ 3 eV 2 


tan 2 e 


0.23 


0.35 


tan 2 6atm 


0.83 


1.0 


tan 2 c hooz 


3.3 x 10~ 2 


<2.6 x 10~ 2 


g.o.f. 


3.38 





Table 1. Best fit to conventional experimental data. All masses are running masses at 1 GeV 
except the top quark mass which is the pole mass. Note that we use the square roots of the 
neutrino data in this Table, as the fitted neutrino mass and mixing parameters (m), in our 
goodness of fit (g.o.f.) definition. 



Family Replicated Calculation of Baryogenesis 



159 



in the decay of one N R i is just ei. The total decay rate at the tree level is given by 



I Nit 



47t(cfj v 



Mi , 



(12) 



where M° can be expressed through the unitary matrix diagonalising the right- 
handed neutrino mass matrix Vr: 



M° = M°Vr , 
V R M R M{ V R =diag ( M?, M 2 , M 2 

The CP violation rateis computed according to 1 8 9 1 



(13) 
(14) 



47t(4> ws ) i ((MD)tMD; 



(15) 



where the function, f (x), comes from the one-loop vertex contribution and the 
other function, g(x), comes from the self -energy contribution. These e's can be 
calculated in perturbation theory only for differences between Majorana neutrino 
masses which are sufficiently large compare to their decay widths, i.e., the mass 
splittings satisfy the condition, |Mt — Mj| 3> |Ii — fj|: 



f(x) 



II +x)ln- 



9W 



1 -x 



We as usual 1 2J introduce the dacay rate relative to 



Ki 



2H 



Mpi anc k 



T=Mi 



1.66(4) 



Mi 



(1 = 1,2,3) 



(16) 



(17) 



where r\ is the width of the flavour i Majorana neutrino, Mi is its mass and i 
is the number of degrees of freedom at the temperature Mi (in our model ~ 1 00). 

In order to estimate the effective K factors we first introduce some normal- 
ized state vectors for the decay products: 



k=l 

([M°(Mi 



ki 



M°(Mi] 



2i 



M°(Mi; 



3i 



Then we may take an approximation for the effective K factors: 
K effl =K 1 (M 1 ) , 

K eff2 =K 2 (M 2 ) + |(2|3)| 2 K 3 (M3) + |(2|1)| 2 ^{Mi 



K 



efl3 



K 3 (M 3 )+|(3|2)| 2 K 2 (M 2 ) + |(3|1)| 2 K, (M, ; 



(18) 
(19) 
(20) 
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5 Result for baryogenesis 

Using the Yukawa couplings - as coming from the VEVs of our seven different 
Higgs fields - the numerical calculation of baryogenesis were performed using 
our random order unity factor method. In order to get baryogenesis in Fukugita- 
Yanagida scheme, we calculated the see-saw neutrino masses, K effi factors and 
CP violation parameters using N = 10,000 random number combinations and 
logarithmic average method: 



Mi 


= 2.1 x 10 5 


GeV 


, K e(f i 


= 31.6 , 


lei 1 


= 4.62 x 10- 


-12 


M 2 


= 8.8 x 10 9 


GeV 


, K e(f2 


= 116.2 


, Ie 2 | 


= 4.00 x 10- 


-6 


M 3 


= 9.9 x 10 9 


GeV 




= 114.7 


, |e 3 | 


= 3.27 x 10- 


6 



The sign of e\ is unpredictable due to the complex random number coefficients 
in mass matrices, therefore we are not in the position to say the sign of e's. Using 
the complex order unity random numbers being given by a Gaussian distribution 
we get after logarithmic averaging using the dilution factors as presented by |2 2| 

Yb =2.59+^x10-" , (21) 

where we have estimated the uncertainty in the natural exponent according to 
the ref. CHI to be 64 % • /fO » 200 %. 

The understanding of how this baryon to entropy prediction Yb comes about 
in the model may be seen from the following (analytical) estimate 

Y ~ 1 X M 3 ~ 1 in-? (77) 

^r^'M^l' 10 (22) 

where we left out for simplicity the In K factor in the denominator of the dilution 
factor k and where M3 is the mass of one of the heavy right-handed neutrinos in 
our model M3 w (4>b-l)- Since the atmospheric mass square difference square 
root y/Am?~ « 0.05 eV » (4> W s) 2 (WT) 2 /M 3 we see that keeping it leaves us 
with the dependence 

<<Kvs) 2 X ^xlO-.^* (23 ) 



3^0.05 eV - g*Mpi anck W 2 T 4 5 yg7W 2 T 4 



6 Problem with wash-out effects? 



To make a better estimate of the wash-out effect we may make use of the calcu- 
lations by [11 J by putting effective values for the see-saw neutrino mass M and 
fa. The most important wash-out is due to "on-shell" formation of right-handed 
neutrinos and only depends on K or the thereto proportional fa, but next there are 
wash-out effects going rather than by K or m as Mm 2 . In the presentation of the 
results by 1 11 [ fixed ratios between right-handed neutrino masses were assumed. 
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However, in reality a very important wash-out comes form the off-shell inverse 
decay and that goes as 

M ^W^ With ^~ M- (24) 
j 1 ' 

Here we use the notation with fn.j from 1111 : fn.j w Kj • 2.2 • 10~ 3 eV. 

Using such a term (see eq. I24t with the ansatz ratios used in 1111 . M| — 
10 6 Mf and M| = 10 3 Mf one gets for eq. O w 10 6 Mi fru|, while we would 
with our mass ratios (eq. E} Mj w 1/4 ■ 10 10 Mf and Mf « 1/4 • 10 10 Mf 
obtain correspondingly 2 • 10 5 GeV • 1/4 • 10 10 fnj » 1/2 • 10 15 GeV fnj, which 
then being identified with 1 6 Mi use mj would lead to that we should effectively 
use for simulating our model the mass of the right handed neutrino - which is 
a parameter in the presentation of the dilution effects in|llj-Mi use = 1/2 • 
10 15 GeV/10 6 = 1/2 • 10 9 GeV. Inserting this Mi use value for our estimate m 2 ~ 
fn.3 w 0.1 eV gives a dilution factor k w 10~ 4 , a factor 500 less than what 
we used with our estimate using the K e ff's. (Our fn.3 = fn.2 are surprisingly large 
compared to the \J Am^ tm because of renormalzation running .) Using the better 
calculation of [11J which has a very steep dependence - a fourth power say - as 
function of m our uncertainty should also be corrected to a factor 100 up or down. 
So then we have one and a half standard deviations of getting too little baryon 
number. 



7 Conclusion 

We calculated the baryon density relative to the entropy density - baryogene- 
sis - from our model order of magnitudewise. This model already fits to quark 
and lepton masses and mixing angles using only six parameters, vacuum expec- 
tation values. We got a result for the baryon number predicting about a factor 
only three less than the fitting to microwave background fluctuations obtained by 
Buchmuller et al. |12[, when we used our crude K e ff's approximation. However, 
using the estimate extracted from the calculations of [11 1 we got three orders of 
magnitude too low prediction of the baryon number. This estimate must though 
be considered a possibly too low estimate because there is one scattering effect 
that is strongly suppressed with our masses but which were included in that cal- 
culation. But even the latter estimate should because of the steep dependence of 
the result on the parameters be considered more uncertain and considering the 
deviation of our prediction only 1 .56o~ is not unreasonable. 

Since we used the Fukugita-Yanagida mechanism of obtaining first a lep- 
ton number excess being converted (successively by Sphalerons) into the baryon 
number, our success in this prediction should be considered not only a victory for 
our model for mass matrices but also for this mechanism. Since our model would 
be hard to combine with supersymmetry - it would loose much of its predictive 
power by having to double the Higgs fields - we should consider it in a non 
SUSY scenario and thus we can without problems take the energy scale to infla- 
tion/reheating to be so high that the plasma had already had time to go roughly 
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to thermal equilibrium before the right-handed neutrinos go out-of-equilibrium 
due to their masses. We namely simply have no problem with getting too many 
gravitinos because gravitinos do not exist at all in our scheme. 

Another "unusual" feature of our model is that the dominant contribution to 
the baryogenesis comes from the heavier right-handed neutrinos. In our model 
it could be arranged without any troubles that the two heaviest right-handed 
neutrinos have masses only deviating by 1 0% namely given by our VEV param- 
eters X- This leads to significant enhancement of the ei and £3 which is crucial 
for the success of our prediction. There is namely a significant wash-our taking 
place, by a factor of the order of k = 10~ 3 to 10~ 6 . It is remarkable that we have 
here worked with a model that order of magnitudewise has with only six ad- 
justable parameters been able to fit all the masses and mixings angles for quarks 
and leptons measured so far, including the Jarlskog CP violation area and most 
importantly and interestingly the baryogenesis in the early Universe. To confirm 
further our model we are in strong need for further data - which is not already 
predicted by the Standard Model, or we would have to improve it to give in prin- 
ciple accurate results rather than only orders of magnitudes. The latter would, 
however, be against the hall mark of our model, which precisely makes use of 
that we can guess that the huge amount of unknown coupling constants in our 
scheme with lots of particles can be counted as being of order unity. 
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Neutrino Oscillations in Vacuum on the Large 
Distance: Influence of the Leptonic CP-phase. 
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ITEP, Moscow, 117259, B.Cheremushkinskaya, 25 

Abstract. Vacuum neutrino oscillations for three generations are considered. The influ- 
ence of the leptonic CP-violating phase (similar to the quarks CP-phase) on neutrino oscil- 
lations is taken into account in the matrix of leptons mixing. The dependence of probabili- 
ties of a transition of one kind neutrino to another kind on three mixing angles and on the 
CP-phase is obtained in a general form. It is pointed that one can reconstruct the value of 
the leptonic CP-phase by measuring probabilities for a transition of one kind neutrino to 
another kind averaging over all oscillations. Also it is noted that the manifestation of the 
CP-phase in deviations of probabilities of forward neutrino transitions from probabilities 
of backward neutrino transitions is an effect practically slipping from an observation. 

1 Introduction. 

It is unclear up to now in spite of great number of papers devoted to the investi- 
gation of neutrino oscillations, what is the real precision of experimental values 
of three mixing angles and masses of neutrino from different generations? And 
consequently do neutrino really oscillate? The central values of these angles and 
values of errors obtained in different papers and given in our paper change from 
author to author and from paper to paper. Therefore these data are very suspi- 
cious. Below in this paper we gave the value of this precision approximately be- 
cause it is defined very roughly. Nevertheless the investigations of neutrino look 
rather encouraging since the set of large perspective devices (K2K in Japan II 121 , 
CERN-GRAND Sacco (CNGS) in Europe and Fermilab-Soudan in USA) and 
some small but also perspective devices in another regions of the Earth began to 
work recently or will begin to work in the near future. In particular, the preci- 
sion of defining of the values of v T and "v^ masses and also the values of sines 
of the neutrino mixing angles will be appreciably improved in the nearest future 
(in one or two years). Furthermore, we believe that attempts to obtain the value 
of the CP-phase from experimental data will be made in the future in spite of 
apparent present-day hopelessness. 

The present work is devoted to neutrino oscillations and, in particular, to a 
possible manifestation of the leptonic CP-phase in neutrino oscillations. At the 
second section of this paper we consider the standard theory of the neutrino 
oscillations with regard for the leptonic CP-violating phase. Then we give the 

* ryzhikh@heron.itep.ru 
** termarti@heron.itep.ru 
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formulas for the probabilities of the conservation of the neutrino kind and for 
the probabilities of the neutrino transition to another neutrino kind with some 
examples of the possible manifestation of the leptonic CP-phase using modern 
experimental data. And then we investigate the difference between the y a — > -vp 
transitions probability and the "Vp — > -v a transitions probability and the possible 
influence of the leptonic CP-phase on this difference. 



2 Standard theory of neutrino oscillations with regard for 
leptonic CP-violating phase. 

In this section we describe the standard theory of neutrino oscillations including 
the leptonic CP-phase. So, neutrino ("Ve)i_, hv)i_, ("Vt)l which were born in the 
decay reactions or in collisions do not have definite masses. They are superposi- 
tions of neutrino states "Vi ,^2,^3 with definite masses, and their wave functions 
are: 

3 

-vp(x,t) = 2jv- l ) pk ^ k (x,t), P = e,n,T, k = 1,2,3. (1) 

k=l 

Here it is supposed that -v^ = (~vi ,^2,^3) are the wave functions of the neutrino 
with definite masses moved in a beam along the axis OX with not small momen- 
tum Ip-yl > irtv and ultrarelativistic energy E k = \Jvl + — Ip-vl + Tn£/2p-v, 
k = 1,2,3. Thus their wave functions look like: 

y k (x, t) = e^'e^'VlO) = e^S^O) (2) 

Mixing of leptons, i.e. mixing of neutrino when the mass matrix of "elec- 
trons" of three-generations is diagonal, is defined by a unitary 3x3 matrix V* 1 = 
^mns Maki-Nakagava-Sakata. This matrix depends on three mixing angles of 
the leptons $12, #13 an d $23- It is similar to CKM matrix of quarks mixing and 
has a well-known form: 

(C12C13 S12C13 si 3 e- i6l \ 

-S12C13 -Cl2S23Sl 3 e i61 C12C23 - Sl 2 S23Sl 3 e i61 S23C13 (3) 
S12S23 -Cl 2 C23Sl 3 e i51 ~Cl 2 S 2 3 ~ Si 2C23SI 3 e" 1 C23C1 3 / 

Note that the V* 1 can be represented in the form of the product of three matrices 
of rotations (or of mixing of two generations in pairs) 612, 613(61) and 623. It is 
easy to verify that V 1 = 612613(60623/ where: 

ci 3 0si 3 e- i5 i 
1 
-si 3 e l61 C13 

(4) 

'1 

O23 = I C 2 3 S 2 3 

-S23 C 23 



/ C12 S12 0\ 




-S12 C12 


,6 13 (6 l )= j 


V 01/ 
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here and in 6; is the leptonic CP-violating phase. Its value is not known up 
to now, sometimes, for example, it is considered to be equal to whereas its ana- 
logue - the quark CP-phase 6 q seems to be close to 7t/2 |4J. Acting by matrix J3J 
'v, 

on column ^ = [ y 2 I we obtain according to (JJ: 




[C12C13-V1 (0) + Si2Ci3V 2 (0)e- i( P^ 

m 2 

[-(S12C23 + Cl 2 S23Sl3e i6l )^1 (0) 

+ (c 12 c 23 -s 12 s 2 3Si3e i6l N2(0)e- i( P 21 
+ c 13 s 23 'V3(0)e- i< P3i]e-^ t 

[(S12S23 - C 12 C23Sl3e i51 )-vi (0) 

- (C12S23 + si2C23Si 3 e i5l h2(0)e- i( P 2 ' 

m 2 

+ ci 3 C23^3(0)e- i( P^]e^^ t 
where, using dependence 10 of neutrino states on the time t = L/c we have: 
imr - mrl . .jm; -m 2 )(eV 2 ) 



(5) 



2p^ 



m?) (m? 
'-t = 1.27- l 



E v (MeV) 



-L(m) 



(6) 



where E-v ~ cp^, is an energy of the neutrino beam: E-v » m.3 > m.2 > mi . 
Neutrino states with definite masses are mutually orthogonal and are normalized 
to unity. Using these statements we can easy obtain expressions for probabilities 
of a transition in vacuum of one kind neutrino to neutrino of another kind during 
the time t. 

For probabilities of conservation of e, p., T-neutrino kind we have, respec- 
tively: 



( P(v e v e 



1-2 r 2 1 <.2 r 2 ptcpn 



5 2 p iip3i [2 



IIC13S12 + C 12 Sl3S23e 

+ cf 3 si 3 e i<P3 'l 2 

P(^t^t) = IIS12S23 - Cl 2 C23Sl3e l6l | 2 

+ cf 3 c2 3 e i ^ 1 |2 



i6l l 2 + |c 12 C23 - s 12 S23Si 3 e i5l | 2 e ic P 2 



IC12S23 + S12C23S13e i6l | 2 e i^p2, 



(7) 
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And for probabilities of transitions of y a neutrino to neutrino of another kind "Vp 
we obtain: 

' P(^e^ H ) = |Cl2Cl3(Cl3Sl2 + Ci 2S23S1 ) 

- Cl3Sl2(Cl2C23 - Si 2S 2 3 Si 3 e lSl )e 1<P2 1 

-Sl3Cl 3 S23e i(5l + Cp3l) | 2 
P(^e^T) = |C12C13(S23S12 - Cl 2C23S1 3e lSl ) 

-Cl 3 Sl 2 (Cl2S23 +C23Sl2Sl3e i • 5l )e ilp2, (8) 

+ Sl3Cl3C2 3 e i(6l + ip ^ ) | 2 
P(^H^t) = I(C13S12 + Cl2Sl3S23e l5l )(Sl2S23 - Ci 2C23S1 3e~ l61 ) 

+ (C12C23 - si2Si 3 S23e i5 )(ci2S23 + C23SI 2Sl 3 e i5 )e i(p2 ' 

-C2 3 cf 3 S23e ilP31 | 2 



3 Probabilities of the change of the neutrino kind 1 — P(^v a ^a) 
and of neutrino ~v a transition to neutrino "Vp of another kind: 

After not complicated, but cumbersome transformations of the formulas 
we have complete expressions for the probability of the change of neutrino kind 
1 — P(^a"v«) and for the probabilities of transition of one kind neutrino to neutrino 
of another kind. But these formulas are very complex for an analysis and because 
of experimental peculiarities of the neutrino registration it is more convenient to 
use probabilities averaging over oscillations, i.e. over phases of neutrino of 
the continuous spectra. Therefore we adduce these complete formulas only for 
references. 

1 -P{y e y e ) =cf 2 sin 2 (2Si3)sin 2 (cp3i/2) + c4 3 sin 2 (2Si 2 ) sin 2 (cp 2 , /2) 
+ s 2 2 sin 2 (2fli 3 )sin 2 (cp32/2) 



1 -Phvvn) ={c43sin 2 (2^ 2 ) + s4 2 s 2 3sin 2 (2^ 2 3) + s43s43sin 2 (2^i 2 ) 

+ c| 2 s 2 3 sin 2 (2-& 23 ) +cos6isin(4{h 2 ) sin(2-923)(si3C 23 - sf 3 s| 3 ) 

- cos 2 6 l sf 3 sin 2 (2^ 23 ) sin 2 (2#i 2 )}sin 2 (cp 2 i /2) 

+ {s 2 2 c 2 3 sin 2 (2^ 2 3 ) + c 2 2 s^ 3 sin 2 (2di 3 ) 

+ cos 6 v s^ 3 ci 3 sin(2^i 2 ) sin(2£ 23 ) sin(2#i 3 ) } sin 2 ( cp 3 1 /2) 

+ {c 2 2 c 2 3 sin 2 (2S 23 ) + s 2 2 s^ sin 2 (2di 3 ) 

-cos6 l s^ 3 Ci3 sin(2-&i 2 )sin(2^ 23 ) sin(2-&i3)}sin 2 (cp 32 /2) 

1 - P(t x t t ) = {s* 3 sin 2 (2di 2 ) + s{ 2 s 2 3 sin 2 (2^ 23 ) 

+ c 4 23 s 4 u sin 2 (2di 2 ) + c? 2 s 2 3 sin 2 (25 23 ) 

+ cos6 l sin(45 12 ) sin(2^ 23 )(s^ 3 c^3 - s 2 3 s^ 3 ) 

- cos 2 6 l s 2 3 sin 2 (2^ 23 ) sin 2 (2^i 2 )}sin 2 (cp 2 i /2) 
+ {s 2 2 c 2 3 sin 2 (25 23 ) + c 2 2 c4 3 sin 2 (2^i 3 ) 

- cos 6tc| 3 ci3 sin(2-&i 2 ) sin(2-& 23 ) sin(2-&i 3 )} sin 2 (cp 31 /2) 
+ {c 2 2 c 2 3 sin 2 (2S 23 ) + s 2 2 c4 3 sin 2 (2^i 3 ) 

+ cos 6tc| 3 ci 3 sin(2-0i 2 ) sin(2-0 23 ) sin(2-&i 3 )} sin 2 ( 932/2) 
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PK^V) = {{sin 2 (25i 3 )(s 2 3 +0^3 + s| 2 s 2 3 ) 

+ jCm sin(25i 3 ) sin(25 23 ) sin(45i 2 ) cos 60 
- 2c 2 3 sin 2 (25 12 )(c| 3 - sf 3 s| 3 ) cos((p 2 i ) 
-2s 2 3 sin 2 (2di 3 )( C 2 2 COS((p3l ) + S 2 2 COS((p3 2 )) 
+ C13 sin(25i 2 ) sin(25i 3 ) sin(25 23 ) 

• (s 2 2 cos(6 v + CP21) -c 2 2 cos(6 v - cp 2 i)) 
+ C13 sin(25] 2 ) sin(25i 3 ) sin(25 23 ) 

• (005(6! + (p 32 ) - cos(6 l - cp 31 )) 
+ 2c 2 3 c 2 3 sin 2 (25i 2 )} 



P(w) = J{sin 2 (25 13 )(c 2 3 + C 4 2 C 2 3 + S 4 2 C 2 3 ) 

- jC U sin(25i 3 ) sin(25 23 ) sin(45i 2 ) cos6i) 
+ 2cf 3 sin 2 (25 12 )(s| 3 - c 2 3 s| 3 ) cos(cp 2 i ) 



-2cf 3 sin 2 (25i 3 )(c 2 2 cos(<p 3 i) + s 2 2 cos(cp 32 )) 
+ C13 sin(25i 2 ) 8111(25! 3 ) sin(25 23 ) 

• (cf 2 C0S(6! - (p 2 1 ) - S 2 2 C0S(6! -+ 

+ ci 3 8111(25! 2 ) sin(25i 3 ) sin(25 23 ) 

• (cos(6! + (p31 ) - cos(6! + cp 32 )) 
+ 2c 2 3 s 2 3 sin 2 (25i 2 )} 



<P2l)) 



P(-v^t) = zf{ 2s is sin 2 (25i 2 )cos 2 (25 23 ) 



+ (C^ 3 + C 

-[2s 2 ,(c 4 



12 



'12 



-12 



+ sj 2 )sj 3 



sin 2 (25 23 ) 



'13^23 



+ s 4 3 )sin z (25i 2 ) + [2s 2 3 (c 12 



'12' 



- (1 + s? 3 ) sin 2 (25, 2 )] sin 2 (25 23 )] cos((p 2 i ) 
-[2c 2 3 (s 2 2 + c 2 2 s 2 3 )sin 2 (25 23 ) 

2C13 sin(25i 2 ) sin(25i 3 ) sin(45 23 ) cos6!] cos(cp 3 i 



[2c 2 3 sin 2 (25 23 )(s 



2 _2 _ r 2 
12 b 13 L 12' 



- 2C13 sin(25i 2 ) sin(25i 3 ) sin(45 23 ) cos 6 J cos(cp 32 ) 
+ 2ci 3 sin(25i 2 ) sin(25i 3 ) sin(25 23 ) 

• sin6!sin(cp 2 i/2) cos( ( P 31 ^ ^ 32 ) 
+ si 3 sin(45i 2 ) sin(45 23 ) cos 6J1 + sf 3 ] sin 2 (cp 2] /2) 

— C13 sin(25i 2 ) sin(25i3) sin(25 2 3) sin6!sin(cp 2 i ) 

- 2s 2 3 sin 2 (25i 2 ) sin 2 (25 23 ) cos(26!) sin 2 (cp 2 i /!)} 

Note that the probability of -v e y e oscillations does not depend on the lep tonic 
CP-phase in contrast to another probabilities of the neutrino oscillations. After 
averaging these formulas over all phases cptj and taking into account (cos((ptj ± 
= 0, (sin 2 cptj ) = (cos 2 cp^) = 1/2 we have: 

( (1 -P(-Ve-Ve)) = A ee 

(1 - Phwn)) = + B^ cos &i + C w cos 2 6i 
(1 — P("v t Vt;)) = A TX + B TT cos 8i + C TT cos 2 5i 

(P('Ve'Vn)) = A eix + B e ^ COS 6 L 

(P(^ e v T )) = A eT + B eT cos 6 V 
(Phvv T )) = + B^ cos 5i + C HT cos(26!) 



(9) 
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A„ 



A, 



1 

,2 
'13 



cf 3 sin 2 (2* 12 )+sin 2 (2* 13 ; 



+ c 



12 



st 2 ) 



'13 



sj 3 sin 2 (2*, 2 ) + sin- 



sin 2 (2* 23 ) 
! (2-& 13 ))s4 3 



+c2 3 sin 2 (2* 12 ) 



C 



Ifr 2 

2 ^ C 23 



2 2 
~~ S 23 S 13 



sin(2* 23 ; 



sf 3 sin 
sin 2 (2*i 2 



Sl3 

sin(4* 12 ) 
(2* 23 ) 



A TT = l[(c 2 3 + (cf 2 + s 4 2 )s 2 3 )sin 2 (2* 23 ) 
+ (s 4 3 sin 2 (2* 12 ) +sin 2 (2* 13 )) C23 



+s 23 sin 2 



13J 

(2*12) 



-jsi3 sin(2* 23 ) 

' ( S 23 — C 23 S 13) 

— -jS 2 3 sin 2 
•sin 2 (2*! 2 ) 



sin(4* 12 ) 

(2*23) 



A e ^ = 1[(1 + cj 2 + s 4 2 )s 2 3 sin 2 (2* 13 )] B evi = l Cl3 sin(2* 13 ) sin(2* 23 ) 
+2c 2 3 c 2 3 sin 2 (2* 12 )] -sin(4* 12 ) 

A eT = I[(1 + c 4 2 + s 4 2 )c 2 3 sin 2 (2* 13 )] B eT = -l Cl3 sin(2* 13 ) 

+2c 2 3 s 2 3 sin 2 (2* 12 )] • sin(2* 23 )sin(4* 12 ) 

A, T = l[2s 2 3 sin 2 (2* 12 )cos 2 (2* 23 ) B »" = + s 2 3 )s 13 sin(4* 12 ) 



. v,, ,,, i |,| sin(4* 23 ) 

■sin (2*23 H(CT 2 + Su) s i3 c - 1 s 2 sin 2 f2*nl 

_Lr 4 + r 4 — 1— q 4 11 HT ~~ ~4 S 13 Sm 1 Z ^12J 

+ C 13 + C 12 + S l 2 iJ „;„2f 



•sin z (2* 23 ) 

(10) 

These expressions are organized in such a way for to emphasize the influence of 
the leptonic CP-phase on the averaging probabilities of the neutrino oscillations. 

Note that the probabilities of the change of neutrino kind and the probabil- 
ities of transitions to another two neutrino states obviously obey the following 
rules: 

1 - P(-v a -v a ) = Ph^-vp) + P(^a^ T ), where a, (3,y = e, 



4 Some examples of the manifestation of the leptonic CP-phase 

In this section we give some examples demonstrating a possible dependence of 
the probabilities on CP-phase. But first, for the convenience we introduce new 
designations: 

bik = Bik/Aik, Cik = Cik/Atk 

The first set of possible values for mixing angles taken from experimental 
data |5| is: 
example a) 

*u = (42 ± 2)°, * 13 = (4.0 ± 0.5)°, * 23 = (43.6 ± 0.5)° (11) 

(small mixing of 1,3 generations was obtained from experimental data of nuclear 
reactor CHOOZ |3 [). Here and below in our paper an average error is taken from 
tables adduced in papers 1 5 6 1 . These values are preliminary and are used below 
for our estimations. In this case all coefficients btk, Ctk in formulas have very 
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small values, in particular because of smallness of mixing angle S13 = sin $13. 
The table of all coefficients of formulas for this case is: 





= 0.499; 






A w 


= 0.636, bjxji = 


= 0.0058, = 


-0.0038; 


( A TT 


= 0.613, b TT = 


0.0055, c TT = - 


-0.0040; 




= 0.261, b e ^ = 


0.014; 




A eT 


= 0.238, b eT = 


-0.015; 






= 0.373, b^ T = 


0.0005, c^ T = 


-0.0032. 



(12) 



As we can see, the ratio of the number of the u-neutrino to the number of the 
T-neutrino produced in the initial beam of electron neutrino "V e at a large distance 
from the source is: 



(PCveVn)) _ A^ (1 +b e , x cos6 l ) 



A 



en 



where b e 



A„ 



(1 +b eT cos60 A e 



1 + (ben ~~ b eT) cos 61 



(13) 



2b 



2.8%, (as b e 



>e.\i 



), with 4^ 



1.1. Thus, the 



contribution of terms containing cos 61 to the ratio of probabilities (P("v e 'V|j.)) an d 
(PtVe^x)) is of order of 3%. Therefore an experimental observation of the CP- 
violating phase manifestation is very difficult for this set of mixing angles (see 
Fig.l). 

The second set of possible values for mixing angles taken from experimental 
data is: 
example b) 



*12 = (42.0 ±2.0)°, d 



13 



(14.0 ± 1.0)°, -9 23 = (43.6 ±0.5)°, 



(14) 



(this example is not in a good agreement with experimental data |3| because, 
although it gives appropriate values for {h 2 and -823 ^llt . the value of fh 3 is rather 
large here |7|). 

In this case, i.e. at sin{>i3 = 0.24, coefficients b^, cue have larger values in 
comparison with the previous case. The values of these coefficients are of the 

8.4%. So, we have in 



2b 



order of several percent, in particular, b e[ x — b f 
this case: 

A ee = 0.548; 

A HH = 0.645, b^ = 0.019, = -0.044; 
A TT = 0.627, b TT = 0.017, c TT = -0.045; 
A e ^ = 0.283, b e ^ = 0.041; 
A eT = 0.265, b eT = -0.043; 



(15) 



A HT = 0.348, b HT = 0.0077, c^ T 



-0.0409. 



Since b e |x — b 



t 8.4%, the contribution of terms containing cos 5; to the ratio 
of probabilities (P("v e "V)) and (P("Ve~v-r)) is of order of 8.4%. The results of corre- 
sponding measurements seem to be very interesting (see Fig.2). 
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m 

Q-634 
0.633 




tt/2 

0.265 
0.264 
0.263 
0.262 
0.261 
0.26 
0.259 
0.258 



tv 3tv/2 2tv 




<1-P(i/ t i/ t )> 

0.614 
0.613 
0.612 
0.611 
0.61 
0.609 
0.608 
0.607 

tt/2 

<^(^r)> 

0.242 
0.241 
0.24 
0.239 
0.238 
0.237 
0.236 
0.235 
0.234 




3tt/2 




<i>(^)> 




Fig. 1. In this figure the dependences of the probabilities averaging over neutrino oscil- 
lations on the leptonic CP-phase are shown. These dependences are small and difference 
between maximal and minimal values of probabilities equals approximately 0.7% in the 
y^yi and v T -v T cases, 0.8% in the -v e "V and v e "v T transitions and 0.03% in the y^y T tran- 
sitions. 



5 The difference between the ~v a — > -v p transitions probability 
and the — > ~v a transitions probability and leptonic 
CP-phase 

Let us consider -v a — > "Vp transitions in neutrino oscillations (denote them "for- 
ward" transitions) and compare them with "Vp — > -v a transitions (denote them 
"backward" transitions). At 6; ^ the probabilities of forward transitions differ 
from probabilities of backward transitions. Note right away that the probability 
of the forward transition coincides with the probability of the backward transition 
after averaging over all oscillations, while before averaging these probabilities are 
different provided CP-violating, i.e. at by ^ 0. And the difference between these 
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0.65 
0.645 
0.64 
0.635 
0.63 
0.625 
0.62 
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0.61 
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0.294 

0.29 
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0.282 

0.278 

0.274 

0.27 




0.362 
0.358 
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0.35 
0.346 
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0.338 
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(1-P(v r u r )) 

0.63 
0.625 
0.62 
0.615 
0.61 
0.605 
0.6 
0.595 
0.59 
0.585 

tv/2 

0.278 
0.274 
0.27 
0.266 
0.262 
0.258 
0.254 



7T 3tt/2 2tt 




tv 3tv/2 2tv 




Fig. 2. In this figure the dependences of the probabilities averaging over neutrino oscilla- 
tions on the leptonic CP-phase are not so small as in the previous case (Fig.l.) because 3 
is not small and difference between maximal and minimal values of probabilities equals 
approximately 4% in the 7^7^ and 7 T 7 T cases, 2.5% in the 7 e 7 M and -v e "VT transitions and 
2.8% in the 7^ transitions. Thus the role of -9] 3 in the manifestation of leptonic CP-phase 
in neutrino oscillations is very important. 



probabilities is proportional to sin 61. These statements are direct consequences 
of general formulas for P('v e 'v lx ), P^e^t) and P( , v ( /v-r) which given above. So, 
on obtaining probabilities of backward transitions by replacement 6; — > —61 we 
subtract probabilities of forward transitions from them. As a result we have: 

Phv^e) - Pt^e^v) =a (sincp 2 i +sincp 32 - sincp 31 ) sinS,. 
P(v T v e ) -P(v e v T ) =-a (sincp2i + sin cp 32 - sin cp 31 ) sin6i (15) 

P("VtV) - P("VVt) = ao(sincp2i -2sin^cos (cp31 + <P32) )sin6 l 
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where ao = 2 c i3 sin2£i2 sin2£i3 sin2-923- Here phases (P21 , CP32, cp3i depend on 
the time of neutrino flight in vacuum t (see 0) or, in other words, on distance 
(base-line) L between points of neutrino birth and neutrino absorption and also 
on difference of squared masses Am?j = m? — m.?. For the most possible values 
of neutrino masses and theirs average errors taken from experiment |5 6[: 

m 3 = (1/17 ± 1/50)eV, m 2 = (1/175 ± 1/300)eV, mi < m 2 (17) 

We can see that the distance of the neutrino oscillations, i.e. the base-line for the 
experimental devices must not be less 1 3 m. The experimental definition of the 
CP-phase based on the correlation H6\ would be the most natural, however now 
it is practically impossible because beams of different types (for example v e and 
"Vyi) of neutrino (that is, obtained in different reactions) but with the same energy 
are required for the experiment. This problem possibly will be solved in the fu- 
ture, but until now all experimental data were obtained only for beams of the 
neutrino with the continuous energy spectra. The cause of this problem consists, 
in particular, in very small cross-sections of neutrino interactions. 

In current experiments we deal only with probabilities of transitions of neu- 
trino with continuous energy spectra in the initial beam, i.e. with all phases J6j 
averaging over oscillations i^l, J10> which also depend on the CP-violating phase 
Si. The main idea of my talk consists in a suggestion to find the value of the CP- 
phase using data of experiments with large base-line and formulas l9l,fi0l. 

Note that coefficient ao defining the value of the effect of t— symmetry vio- 
lating J16I is not too small in both a) — and b)— cases: 

a) a = 0.07 b) a = 0.23 

Moreover in b-case it is large. Therefore measurements of this effect are possible 
although they are difficult. 

6 Conclusions. 

So, the main results of our work are the following: 

• The expressions for the probabilities of neutrino oscillations were obtained in 
the explicit form with regard for the leptonic CP-phase. 

• The manifestation of the leptonic CP-phase in neutrino oscillations was inves- 
tigated by the example of the probabilities averaging over oscillations. Using 
modern experimental data the model calculations and numerical estimates 
were done. 

• The question of the t— symmetry violation for neutrino oscillation was ana- 
lyzed. And there was established that the difference between the "forward" 
probabilities and the "backward" probabilities was proportional to sine of the 
leptonic CP-phase. 
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Abstract. We show that neutral fermions may have an additional source of time rever- 
sal violation by associating time reversal with gauge group representations of fermions 
through the method of group extensions as in the case of parity. This provides a new source 
of time reversal violation for neutral particles. 



1 Introduction 



Topic of time reversal violation in neutrino processes (especially in oscillations) 
has gained much theoretical attention 1 1 1 and there are proposals to study these 
effects experimentally |2|. The main interest in these studies is to use T- violation 
as the measure of CP violation through fermion mixing matrices. In this study we 
show that the only source of time reversal violation in neutral fermion (hence in 
neutrino) processes is not the one which is related (through CPT theorem) to CP 
violation from complex phases in fermion mass matrices. There is an additional 
possible source of time reversal violation for neutral fermions arising from an- 
other kind of fermion mixing other than the Dirac and Majarona ones as we shall 
see below. 



2 Four-component fermions and parity 

Let us consider two 2-component fermions; xi / X2 which transform under the 
same representation of SL(2,C) 

Xi -fre^-^xi, X2 -> eH -^X2. (1) 

Assume that Xi / X2 belong to different gauge groups or to different representa- 
tions of the same gauge group. The simplest Lagrangian which can be constructed 
out of xi / X2 with related kinetic terms and quadratic interaction terms is 

£ = iX^auD[ 1 1) Xi +ix 2 ^D[ l 2) X2 + m(xUff2X2+xUff2X 1 )+b.c. (2) 

where h.c. stands for taking the complex conjugate of the terms preceeding it and 
oyD u = D a ■ D, D^ 1 1 = 9 U + i 9l B^ 1 >, D^ 2 ' = 9 U + ig z B^ ] with B^ 1 \ being 
* E-mail:erdem@likya.iyte. edu.tr 
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the gauge fields coupling to xi and X2, respectively. We assume an unbroken 
electromagnetic gauge group U(l ) em . So xi / X2 have opposite electric charges. 
Instead of identifying Eq.10 as a Lagrangian of two (two-component) fermions 
interacting through a quadratic interaction term one may pass to a 4-component 
formulation so that Eq.J5} becomes 



C = -tv!>0 m P L ^ + ii];0 (2) P R ^ + Tm|n|> + h.c. (3) 



Here 



* = fe) = fe) ' * =Vt °' Pl = I (1+Ys) ' p * = 2^-75) (4) 



where in the chiral representation employed 

»-(;•) 

The difference between Eq.J5J and Eq.lO is that, in Eq.lO there is a single 4- 
component massive fermion while in Eq.10 there are two interacting 2-component 
massless fermions. So the 4-component formulation simplifies the field theoretic 
calculations but the price to be paid is that the upper and lower two compo- 
nents transform under different gauge interactions. In fact this is one of the basic 
observations behind the parity violation in the standard model of electroweak 
interactions. Mathematically speaking the procedure given above where space 
reflection is associated with a Z 2 group corresponding to the exchange of the 
gauge group representations of left-handed and right-handed fermions is known 
as the method of group extensions |3|. In particular here the procedure corre- 
sponds to the pullback of two Z 2 groups, one associated with space reflection 
and the other one associated with the exchange of the gauge group representa- 
tions of the left-handed and right-handed fermions via the isomorphism relating 
the two Z 2 groups |4|. In the next section we shall employ a similar procedure for 
time reversal. 

A field theoretic formulation of these observations will make the picture 
clearer and more precise. First we take two 2-component fermion fields, xi and 
X2 (coupling through a quadratic interaction), with 1 5 1 



Xi 



X2 



d " p ! [u 1 (p)b 1 (p)e- il " x +vi(p)dUp)e il " x i 



u 2 (p)b 2 (p)e- ip " c + v 2 (p)d 2 (p)e ip ' x ] (6) 



where ui , u 2 , vi , v 2 are 2-component spinors with vi - io~ 2 u';, v 2 - icr 2 U2, u 2 - 
ia 2 u|, v 2 - iff 2 u 2 and bi, b 2 , d^, d\ are the annihilation, creation operators for 
Xi /X 2 an d p • x = E • t — p x. In the 4-component formulation xi /X2 are replaced 
by a single 4-component fermion with 



t|> = 



( f. V 3/2 g= [u(p, (T)b(p, o-)e-^' x + v(p, ff)dt(p, o-Je^T (7) 
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where cr stands for either of helicity or spin. In the chiral basis I5I6I7I 
u(p ) 1/2) = ( U 1 ), u(p ) - 1 /2) = (. cr ^ 

V ^ 1 / 2) K V 1 )' V(P '- V2) = (ia 2 ul) (8 > 
b(p,1/2) =b 1 (p), b(p,-l/2) =d 2 (p), 
d(p,l/2) = di(p), d(p,-1/2) =b 2 (p) 

where ±1/2 stands for different helicity states. In the mass basis (which is the 
solution basis for the Dirac equation) 1 8 [ one may write 




e -tp-x (9) 



[u(p, cT)b(p, a)e- ilp x + v( V , &)tf[p, a^v* 

(Ui +iff2U|)bl + (v 2 + 10-2^2 )dz 

(iff 2 u| - ui )bi 
u 2 -iff 2 u^)b 2 \ gip . x 

which is a sum of two independent spin states. 

The lower and upper two components of 4> in the chiral representation, xl 
and xr are related by space reflection and interchange of xl and xr amounts to a 
space reflection operation plus interchange of the gauge group representations of 
Xl and xr |9|. This is the case because in the 4-component formulation xl and xr 
form a single system ij> so their interchanges are mutually related to each other. In 
the next section we shall investigate such a possibility for time reversal. We shall 
couple two 4-component fermions (with different gauge group representations) 
through a quadratic coupling where each 4-component fermion in the coupling 
can be redefined so that each is the time reversal of the other. Then the result is a 
single system i.e. an 8-component fermion consisting of the original 4-component 
fermions. In this way time reversal is associated with the interchange of the gauge 
group representations of the original two 4-component fermions. This, in turn, 
enables one to break the time reversal invariance in a way other than the complex 
phases in the fermion mixing matrices. 



3 8-component fermions and time reversal 

In the previous section we have seen that one can interprete a quadratic interac- 
tion term between two 2-component massless fermions as a fermion mass term 
by passing to 4-component fermion formulation. In this section we raise the ques- 
tion if one can find a Lorentz invariant quadratic term, other than the mass term 
for 4-component fermions, between two 4-component fermions, which can be in- 
terpreted as a mass term for 8-component fermions. We shall see that this possible 
provided the fermion is neutral and massless ( or almost massless) and the upper 
and lower two components of these 8-component fermions are related by time re- 
versal. This observation opens up the possibility of an additional source of time 
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reversal violation for neutral fermions in a way similar to the case of parity for 
4-component formulation. 

We consider two 4-component fermions given by 

ui , U-2, Vi , V2 are, this time, 4-component spinors with the spinor parts of i|>i , \b 2 
being each others time reversals i.e. 

u 2 (p,ct) ~iy 1 7 3 u*(p,o-), v 2 (p, a) ~ iyVMtP. a) (11) 

The only possible Lorentz invariant quadratic interaction term between 4> 1 and 
ib 2 is 1 10| 

M(^ 2 +^i) (12) 
The simplest possible Lagrangian reads 

C = i-4>T (1 '-ib! -iiI>2Y°^h + rn.iJ;iiI>i + nv\> 2 $2 + MCi|>ty 2 + 4> 2 ^i ) + hx - ( 13 ) 

where$ (1) = y^(9 H -igi B^ 1 '), (2) = Y H ^-ig 2 B[ l 2) ). We take both *h and ^2 
be massless (for example through a chiral symmetry). Then dl3l can be put into a 
simple 8-component form 

C = 'p u ¥ + xW (2) Pd^ + MW + fix. (14) 

where 

V = ( T »0 ) ' Pu = (0 0) ' Pd = (0 i",) 1 

*-(£)• * = * ,r ' r = (uo) (15 > 

So acts as a single fermion whose upper and lower 4-components couple to 
different gauge groups, and the spinor parts of the lower and the upper compo- 
nents are related by time reversal. Because the spinor parts of "ib 1 , \\>2 are related 
by time reversal i|>2 has the same electric charge as i|>i . However because time 
reversal is violated (for example through weak interactions) there will be an ef- 
fective charge asymmetry induced between xi and X2 which forbids i^i^ type 
of terms. This requires that W should be a neutral particle in order to allow ibjij^ 
type of terms (provided the electric charge is a conserved quantity). In summary 
the Lagrangian M4t is applicable only for neutral particles and in that case there 
is an additional source of time reversal violation other than the one due to the 
complex phases in the fermion mixing matrix. If one assigns one of the i|) 1 , ^2 to 
the singlet and the other to a nontrivial representation of the gauge group then 



f27T)3/2 




u 1 (p)b 1 (p)e- i P' x +v 1 (p)dUp)e 



d J p 



(27TP/2 




u 2 (p)b 2 (p)e~ 



ipoc +v 2 (p)d 



P)e l 



(10) 
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the time reversal becomes maximal in the same way as in the standard model of 
electroweak interactions. The propagator of W corresponding to Eq. tl4t is 

where r " = LoL ~ r y) (i6) 



If we take mi = m.2 = m ~ then the mass terms of ipi , \\>z behave as quadratic 
interaction terms which modify the propagator self energy contribution which 
changes Mis to Mis + Y°Tn.r where Ig stands for 8x8 unit matrix. In this case W 
can be still considered as a single fermionic system with some perturbative self 
interaction term. 

The parity and time reversal operators in this case can be infered more eas- 
ily by considering the Dirac and the Dirac-like equations corresponding to the 
Lagrangians © and I14t 

itZ) m P L rl> + mP R \f> =0, i0 (2) P R Tj) + mP L ^ = (17) 
i0 m P u ^ + MP D iJ) =0, i0 U) P D \l) + MP u ^ =0 (18) 

One can determine the parity and time reversal operators by letting x — > — x 
and t — > — t, respectively in the equations given above. One notices that there 
are more than one operator acting on the fields corresponding to each of these 
transformations. For example, for parity one may take either of 

V o\ r0 _fo-i 

for time reversal either of 



vyV K \ p f° 1 
tyV V' V 10 



(20) 



where K stands for complex conjugation of c-numbers on its right hand side. Here 
the first operators in <19t and 1201 correspond to the direct extension of the usual 
parity and time reversal operators for the 8-component case. One notices that the 
second time reversal operator in <2(Jt is unitary At first sight this seems to con- 
tradict with the Wigner's requirement that time reversal is a symmetry of the 
physical systems so it should be anti-unitary |11|. However Wigner's argument 
assumes that the time reversal of an operator corresponding to a physical observ- 
able Q at time t, Q(t) = exp(iHt)Qexp(— iHt) is given by exp(— iHt)Qexp(iHt). 
In other words he implicitly assumes that the Hamiltonian H is invariant under 
time reversal. Because the additional time reversal operator in J20I only applies 
to neutrinos (whose time reversal properties are not well known) and it does not 
leave the Hamiltonian invariant in general (its violation may be even maximal) 
so the Wigner's argument does not apply to the present case. 

While both operators in l !19t or i20\ correspond to the same space-time re- 
flection the physical operations to perform these transformations are different. In 
other words the experiments designed to study parity or time reversal in a phys- 
ical process are different in each case. The parity transformation corresponding 
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to the first operator in J19t is accomplished by interchange of the left-handed and 
right-handed components of the 4-component fermions. Because this interchange 
is associated with space reflection and each of left-handed and right-handed com- 
ponents are associated with opposite helicities one can check parity invariance 
under this operation by checking the space reflection in a given fermionic process. 
However the second operation in J19I does not correspond to a simple exchange 
of the left-handed and right-handed components so it corresponds to preparing 
a different experiment where the state vectors for the fermions are transformed 
accordingly. In other words it is not possible to test the invariance of parity due to 
the second transformation in H9\ by checking space reflection invariance in the 
same fermionic process, i.e. in this case one should first study the original process 
and then the process corresponding to the parity transformed one. A similar ar- 
gument is true for the time reversal operators in Eq. J20> . The first (the one which 
corresponds to the usual time reversal) corresponds to preparing a different pro- 
cess which corresponds to the time reversed of the original process while the test 
of time reversal for the second operator can be achieved, in principle, in the same 
process by comparing the intensity of the events in the direction forward in time 
and the events in the direction backward in time. While this is not applicable in 
practice one may test the time reversal in the same process in this case in indirect 
ways. The spinor part of W, U can be written in the following way 

U=f^y U, = n,:l. U;=("M. (2D 





v u 2 

In the chiral basis the components of U may be related as 1911 

(1) 1 i , 1 (1) (2) 1 , > (2) 

u '= — (po + o-- p)u 2 u =— po-ff-pUj 

m m 

u' 1 ' =-(- Po + o-.p)ui 2 \ 4 1 ' = l(-p - ff .p) u f (22) 
m m 

We notice that both ib i and \p2 have exactly the same helicities. So if two kinds 
of neutrinos are produced which have the same helicities and which do not inter- 
act with one another through gauge interactions (due to Lorentz invariance) one 
may suspect if these additional neutrinos are ip2's (provided xj>i 's are identified 
by the usual neutrinos). If the interaction of Ws is time reversal invariant under 
electroweak interactions then the presence of additional neutrinos (i.e. ij>2's) are 
seen as an extra factor of 2 in the production cross sections of neutrinos when 
compared with their interaction with Z bosons. If there is a time reversal viola- 
tion then this may either be due to the usual source of time reversal violation 
arising from complex phase(s) in fermion mixing matrices or it may be due to 
non-invariance under the interchanges of ipi and \\>2 - Seeking for any asymmetry 
between the gauge interactions of ij>i 's and \\>2's in the same process amounts to 
seeking a time reversal violation due to the second operator in Eq <20> . However 
it may not be easy to disentagle these two types of time reversal violation. 
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4 Conclusion 

We have seen that as in the case of parity one may extend time reversal operator 
by a 7-2 group associated with the interchange of the gauge group representations 
of two coupled fermions whose spinor parts are time reversals of each other. In 
this way one gets the possibility of an additional source of time reversal violation 
for neutral fermions (e.g for neutrinos). This time reversal violation is similiar to 
the parity violation and it may be even maximal as in the case of parity. If this ad- 
ditional time reversal has a comparable magnitude as the usual time reversal vio- 
lation resulting from complex phase(s) in fermion mass matrices then the mixing 
phase can not be direcly derived from the magnitude of time reversal violation 
and vica versa. The presence of the extra time reversal violation can be infered by 
comparing the magnitudes of the time reversal and the usual CP violation in a 
neutrino process (for example in neutrino oscillations) in future experiments. 
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1 Introduction 

I reviewed the general problem of the quark-lepton mass spectrum at the first 
Bled workshop on "What comes beyond the Standard Model" 1 1 1. So, in this talk, 

1 will mainly concentrate on two topics: the Lightest Flavour Mass Generation 
model and the Neutrino Mass and Mixing problem in the Anti-Grand Unification 
Theory (AGUT). 

2 Lightest Flavour Mass Generation Model 

A commonly accepted framework for discussing the flavour problem is based on 
the picture that, in the absence of flavour mixing, only the particles belonging 
to the third generation t, b and t have non-zero masses. All other masses and 
the mixing angles then appear as a result of the tree-level mixings of families, 
related to some underlying family symmetry breaking. Recently, a new mecha- 
nism of flavour mixing, which we call Lightest Family Mass Generation (LFMG), 
was proposed |2[. According to LFMG the whole of flavour mixing for quarks 
is basically determined by the mechanism responsible for generating the physi- 
cal masses of the up and down quarks, m u and ma respectively. So, in the chiral 
symmetry limit, when m u and m.d vanish, all the quark mixing angles vanish. 
Therefore, the masses (more precisely any of the diagonal elements of the quark 
and charged lepton mass matrices) of the second and third families are practi- 
cally the same in the gauge (unrotated) and physical bases. The proposed flavour 
mixing mechanism, driven solely by the generation of the lightest family mass, 
could actually be realized in two generic ways. 

The first basic alternative (I) is when the lightest family mass (m u or rtia) 
appears as a result of the complex flavour mixing of all three families. It "runs 
along the main diagonal" of the corresponding 3x3 mass matrix M, from the 
basic dominant element M33 to the element M22 (via a rotation in the 2-3 sub- 
block of M) and then to the primordially texture zero element Mj 1 (via a rotation 
in the 1-2 sub-block). The direct flavour mixing of the first and third quark and 
lepton families is supposed to be absent or negligibly small in M. 

The second alternative (II), on the contrary, presupposes direct flavour mix- 
ing of just the first and third families. There is no involvement of the second 
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family in the mixing. In this case, the lightest mass appears in the primordially 
texture zero Mn element "walking round the corner" (via a rotation in the 1- 
3 sub-block of the mass matrix M). Certainly, this second version of the LFMG 
mechanism cannot be used for both the up and the down quark families simulta- 
neously, since mixing with the second family members is a basic part of the CKM 
quark mixing phenomenology (Cabibbo mixing, non-zero V C b element, CP vio- 
lation). However, the alternative II could work for the up quark family provided 
that the down quarks follow the alternative I. 

Here we will just consider the latter scenario. 




2.1 Quark Sector 

We propose that the mass matrix for the down quarks (D = d, s, b) is Hermitian 
with three texture zeros of the following alternative I form: 



(1) 



It is, of course, necessary to assume some hierarchy between the elements, which 
we take to be: Bd ~> Ad ~ [bo I S> |<JdI- The zero in the (Md)^ element cor- 
responds to the commonly accepted conjecture that the lightest family masses 
appear as a direct result of flavour mixings. The zero in (Md) 13 means that 
only minimal "nearest neighbour" interactions occur, giving a tridiagonal matrix 
structure. 

Now our main hypothesis, that the second and third family diagonal mass 
matrix elements are practically the same in the gauge and physical quark-lepton 
bases, means that : 

B D = m b + 6 D A D = m s + 5' D (2) 

The components 6d and 6' D are supposed to be much less than the masses of the 
particles in the next lightest family, meaning: 

|5 D |<m s |5' D | <gc m d (3) 

Since the trace and determinant of the Hermitian matrix Md gives the sum and 
product of its eigenvalues, it follows that 

5 D ^ -m d (4) 

while 6' D is vanishingly small and can be neglected in further considerations. 

It may easily be shown that our hypothesis and related equations I0-|4) are 
entirely equivalent to the condition that the diagonal elements (Ad, Bd) of Md 
are proportional to the modulus square of the off-diagonal elements (a D , b D ): 

A D Qd 2 

DD Dd 
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Using the conservation of the trace, determinant and sum of principal mi- 
nors of the Hermitian matrices M D under unitary transformations, we are led 
to a complete determination of the moduli of all their elements, which can be 
expressed to high accuracy as follows: 



^/m-diris 
|M D | = ( ^m d m s m s y/m d m b ] (6) 
^/TridTTib m b — m d y 

Now the Hermitian mass matrix for the up quarks is taken to be of the fol- 
lowing alternative II form: 

(0 c u \ 
M u = A u (7) 
Vu Bu/ 

The moduli of all the elements of M u can also be readily determined in terms of 
the physical masses as follows: 

(0 ^/m u ra t \ 
m c (8) 
^/TUuTat m t — rrtu/ 

The CKM quark mixing matrix elements can now be readily calculated by 
diagonalising the mass matrices M D and M u . They are given by the following 
simple and compact formulae in terms of quark mass ratios: 



l y usl = = 0.222 ±0.004 |V us | exp =0.221 ±0.003 (9) 



rru 



|V cb | = = 0.038 ±0.004 |V cb | exp = 0.039 ±0.003 (10) 



TUb 



|Vubl = = 0.0036 ± 0.0006 |V ub | exp = 0.0036 ± 0.0006 (11) 

As can be seen, they are in impressive agreement with the experimental values. 



2.2 Lepton Sector 

The MNS lepton mixing matrix is defined analogously to the CKM quark mixing 
matrix: 

U = U v U E t (12) 

Here U E and Uv are the unitary matrices which diagonalise the charged lepton 
mass matrix M E and the effective neutrino mass matrix M-y respectively. Assum- 
ing the charged lepton masses follow alternative I, like the down quarks, the 
LFMG model predicts the charged lepton mixing angles in the matrix Ue to be: 
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These small charged lepton mixing angles will not markedly effect atmospheric 
neutrino oscillations, which appear to require maximal mixing sin 20 Qtm ~ 1. 
Similarly, in the case of the large mixing angle (LMA) MSW solution of the solar 
neutrino problem, they are essentially negligible. It follows then that the large 
neutrino mixings should mainly come from the LLv matrix associated with the 
neutrino mass matrix. 

According to the "see-saw" mechanism, the effective mass-matrix M 7 for 
physical neutrinos has the form 

M.-V = -MJM n ' n M n (14) 

where Mm is their Dirac mass matrix, while Mm n is the Majorana mass matrix of 
their right-handed components. Matsuda et al 1 3 [ have extended the alternative I 
LFMG texture to the Dirac Mm and Majorana Mnn matrices. 

The eigenvalues of the neutrino Dirac mass matrix Mm are taken to have a 
hierarchy similar to that for the charged leptons (and down quarks) 

M N3 :M N2 :M N1 ~1 :y 2 :y 4 , y » 0.1 (15) 

and the eigenvalues of the Majorana mass matrix Mm n are taken to have a stronger 
hierarchy 

M NN3 : M NN2 : M NN1 ~ 1 : y 4 : y 6 (16) 

One then readily determines the general LFMG matrices Mm and Mn n to be of 
the type 

/ ay 3 

M N ~ M N3 ay 3 y 2 ay 2 | (17) 
\0 ay 2 1 

and 

/0 |3y 5 

M NN ~M NN3 |3y 5 y 4 fiy 3 | . (18) 
\0 |3y 3 1 

We further take an extra condition of the type 

|A-1|<y 2 (A = a,|3) (19) 

for both the order-one parameters a and (3 contained in the matrices Mm and 
Mnn, according to which they are supposed to be equal to unity with a few 
percent accuracy. Substitution in the seesaw formula M4t generates an effective 
physical neutrino mass matrix M^, of the form: 

IVlv---^ |yll (y -y')-' 1 - (y -y-l I (20) 
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The physical neutrino masses are then given by: 



The predicted values of the neutrino oscillation parameters 




(21) 



sin 2 29 



atm — 



1 , sin 2 29 



2 



3' 



U e3 



2V2 



V, 



Am 2 
Am 2 



atm 



sun 




(22) 



in agreement with atmospheric and LMA-MSW solar neutrino oscillation data. 

The proportionality condition 0, which leads to the LFMG texture, is not 
so easy to generate from an underlying symmetry beyond the Standard Model. 
However Jon Chkareuli, Holger Nielsen and myself have recently shown 1 4 1 that 
it is possible to give a natural realisation of the LFMG texture in a local chiral 
SU(3) family symmetry model. 

3 Fermion Masses in the AGUT Model 

The AGUT model is based on a non-simple extension of the Standard Model (SM) 
with three copies of the SM gauge group — one for each family — and, in the ab- 
sence of right-handed neutrinos, one extra abelian factor: G = SMG 3 x U(l)f, 
where SMG = SU(3) x SU(2) x U(l). This AGUT gauge group is broken down 
by four Higgs fields S, W, T and £, to the usual SM gauge group, identified as the 
diagonal subgroup of SMG 3 . The Higgs field S has a vacuum expectation value 
(VEV) taken to be unity in fundamental (Planck) mass units, while W, T and f, 
have VEVs an order of magnitude smaller. So the pure SM is essentially valid, 
without supersymmetry, up to energies close to the Planck scale. The AGUT 
gauge group SMG 3 x U(l )f only becomes effective near the Planck scale, where 
the i'th proto-family couples to just the i'th SMG factor and U(l)f. The U(l)f 
charges assigned to the quarks and leptons are determined, by anomaly cancel- 
lation constraints, to be zero for the first family and all left-handed fermions, and 
for the remaining right-handed states to be as follows: 



Q f (T T ) = Q f (b R ) = Q f (c R ) = 1 Qf(Hp) = Qf(6 R ) = Q f (t R ) = -1 (23) 



I refer to the review of the AGUT model by Holger and myself at the first Bled 
workshop 1 5 1 for more details. 

The quarks and leptons are mass protected by the approximately conserved 
AGUT chiral gauge charges |6|. The quantum numbers of the Weinberg-Salam 
Higgs field 4>ws are chosen so that the t quark mass is not suppressed, whereas 
the b quark and t lepton are suppressed. This is done by taking the four abelian 
charges, expressed as a charge vector Q = (yi /2, "1/2/2,1)3/2, Qf), for cf^ws to be 
given by: 



Q* ws = Qc R - Qt L = (0,2/3,0,1) -(0,0,1/6,0) = (0,2/3,-1/6,1) (24) 
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We assume that, like the quark and lepton fields, the Higgs fields belong to sin- 
glet or fundamental representations of all the non-abelian groups. Then, by im- 
posing the usual SM charge quantisation rule for each of the SMG factors, the 
non-abelian representations are determined from the weak hypercharge quantum 
numbers yi. The abelian quantum numbers of the other Higgs fields are chosen 
as follows: 

Qw = (0,-1/2,1/2,-4/3) Q T = (0,-1/6,1/6,-2/3) (25) 
Q t = (1/6,-1/6,0,0) Q s = (1/6,-1/6,0,-1) (26) 

Since we have < S >= 1 in Planck units, the Higgs field S does not suppress the 
fermion masses and the quantum numbers of the other Higgs fields W, T, E, and 
$ws given above are only determined modulo those of S. 

The effective SM Yukawa coupling matrices in this AGUT model can now be 
calculated in terms of the VEVs of the fields W, T and £, in Planck units — up to 
"random" complex order unity factors multiplying all the matrix elements — for 
the quarks: 

/WT 2 £, 2 WT 2 £, W 2 T£,\ 
Y u ~ WT 2 £, 3 WT 2 W 2 T (27) 

V £, 3 1 WT J 

/W(T 2 £, 2 WT 2 £,T 3 £A 
Y D - WT 2 f, WT 2 T 3 (28) 

V W 2 T 4 £, W 2 T 4 WT / 

and the charged leptons: 

/WT 2 £, 2 WT 2 £, 3 WT 4 £, \ 
Y E - WT 2 £, 5 WT 2 WT 4 £, 2 (29) 
\WT 5 £, 3 W 2 T 4 WT / 

A good order of magnitude fit is then obtained 1 5 [ to the charged fermion masses 
with the following values for the Higgs field VEVs in Planck units: 

W = 0.179, T = 0.071 £, = 0.099. (30) 

We now consider the neutrino mass matrix in the AGUT model. 



4 Neutrino Mass and Mixing Problem 

Without introducing new physics below the AGUT scale, the effective light neu- 
trino mass matrix is generated by tree level diagrams involving the exchange 
of two Weinberg-Salam Higgs tadpoles and the appropriate combination of W, T, 
£, and S Higgs tadpoles. In this way we obtain: 

, \2 /W 2 £, 4 T 4 W 2 £T 4 W 2 £, 3 T\ 
Mv ~ {<Pws) W 2 £T 4 WT 5 W 2 T , (31) 
Mn \W 2 ^T W 2 T W 2 T 2 ^ 2 / 
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The off-diagonal element (My) 23 = (M.-v)32 dominates the matrix, giving large 
v T — v^ mixing with the following two neutrino masses and mixing angle: 



Although the large mixing angle sin 2 29 RT is suitable for atmospheric neutrino 
oscillations, there are two problems associated with the neutrino masses. Firstly 
the ratio of neutrino mass squared differences Am^/Am^ ~ 2TE, 2 ~ 1 .4 x 10~ 3 , 
whereas the small mixing angle (SMA) MSW solution to the solar neutrino prob- 
lem requires Am^/Am^ ~ 10~ 2 . Secondly the predicted overall absolute mass 
scale for the neutrinos (4>ws) 2 /Mpianck - 3 x 10~ 6 eV is far too small. 

We conclude it is necessary to introduce a new mass scale into the AGUT 
model. Two ways have been suggested of obtaining realistic neutrino masses and 
mixings in the AGUT model: 

1. By extending the AGUT Higgs spectrum to include a weak isotriplet Higgs 
field A with SM weak hypercharge y/2 = —1 and a VEV (A ) ~ 1 eV; also a 
new Higgs field tJ> giving large |J. — T mixing in the charged lepton Yukawa 
coupling matrix is required. 

2. By including right-handed neutrinos and extending the AGUT gauge group 
to G 

extended — (SMG x U(l )b— l)^/ also two new Higgs fields 4>b— l and x 
are introduced to provide a see-saw mass scale and structure to the Majorana 
right-handed neutrino mass matrix. 

Yasutaka Takanishi reported on the second approach |7| at the workshop; so 
I will report on the first approach |8| here. We must therefore consider the intro- 
duction of a new Higgs field i]>, which can yield large mixing from the charged 
lepton mass matrix without adversely affecting the quark mass matrices. With 
the following choice of charges for the i]> field 



m.2 ~ m.3 ~ 



(4>ws) 2 



W 2 T 



sin 2 29 UT ~ 1 



(32) 



Qtp =3Q £ + Qw + 4Q T 



( 



2' 3'6" 




(33) 



we obtain new expresssions for the quark Yukawa matrices: 




WT 2 £, 2 WT 2 £, W 2 T£, 
WT 2 £, 3 WT 2 W 2 T 



(34) 



£, 3 1 WT 




WT 2 f, 2 WT 2 f. T 3 f, 
WT 2 £, WT 2 T 3 
W 2 Ti|; W 2 TEj\> WT 



) 



(35) 




WT 2 f. 2 W 2 T 2 ^ f, 4 ^ 
W 4 T£,tJ) 2 WT 2 
W 3 £, 2 tJ) W 2 T£,il) WT 




(36) 
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As we can see from the charged lepton matrix we will indeed have large mixing 
if (ij>) = 0(0.1 ), so that (Y^n ~~ (Ye)33- In the following discussion we shall take 
4> to have a vacuum expectation value of (t|>) = 0.1 for definiteness. The effect 
of the field \\> on the charged fermion masses is then small, since the elements 
involving ij) do not make any significant contribution to the determinant, or the 
sum of the minors, of the mass matrices, or the trace of the squares YY^ of the 
Yukawa matrices. The mixings of the quarks is essentially unaffected by the terms 
involving \\>, and the only significant effect is on the mixing matrix Ue, which is 
now given by: 




1 0.021 6.4 x 10- 
-0.016 0.75 0.66 
0.014 -0.66 0.75 



where X = ^1 + W 2 J 2 /E, 2 ^ 2 ~ 1 .51 This gives the large mixing required, 

Sin 2 29 a tm - 1, 



(37) 



(38) 



for the atmospheric neutrino oscillations. 

We can further obtain a solution with vacuum oscillations for the solar neu- 
trinos by choosing appropriate charges for the isotriplet Higgs field A. We require 
a large off-diagonal (1,2) element for the neutrino mass matrix and hence we 
choose the charges on A to be 



Qa 



2' '2'3 J 



We then obtain the neutrino mass matrix, 



/ WE, 6 WE, 3 JE, 2 ^ ' 
M v - (A ) WE, 3 W T£iJ) 
\TE, 2 ^ T£iJ) T 2 £,^ 



(39) 



(40) 



This has the eigenvalues, 

(A 



mi 



T£, 2 xh 



T 2 MA 
2 ' 



m 2 ~(A°)W, m 3 ~(A°)(T^ 



T 2 ^ 



(41) 

where the splitting between mi and m.2 comes from the mass matrix element 
(M.-v)33. The neutrino mixing matrix is then given by, 



U, 




(42) 



Hence, using Ue from eqn.|37| we have the lepton mixing matrix U = U^U^,: 



U 



_J_f1 _ X) wt Lfi + X 



WT(1- 



WT(1- 



0.83 -0.016 0.58 
0.38 0.75 -0.55 
-0.43 0.66 0.62 



(43) 
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which, as we can see, has large electron neutrino mixing, as we require for a 
vacuum oscillation solution to the solar neutrino problem. 
We also have the mass hierarchy, 

Ara 2 3 TW 3,^-7 m 
Am 2 3 1 W2 3X10 • (44) 

Hence, if we then take (A ) - 0.2 eV, so that we have an overall mass scale suitable 
for the atmospheric neutrino problem, then we will also have, 

Am 2 3 - 2(A°) 2 T 3 £, 3 xL> 2 - 3 x 1(T 10 eV 2 . (45) 

With such a hierarchy of Am 2 s we effectively have two-neutrino oscillations for 
the solar neutrinos, with the mixing angle given by, 

sin 2 29 sun = 4U 2 ! U 2 3 ~ 0.9. (46) 

So, we have the 'just-so' vacuum oscillation solution to the solar neutrino problem 
with large electron neutrino mixing. We remark that \l e i — —0.016 satisfies the 
CHOOZ electron neutrino survival probability bound (U e 2 is the relevant mixing 
matrix element, since Am 2 2 ~ Am^ S> Am 2 3 ). 

It is also possible to obtain a small mixing angle SMA-MSW solution to the 
solar neutrino problem, with a different choice of charges for A: 

Q A = (-1,4-1,0) (47) 
which gives the quasi-diagonal neutrino mass matrix 

/W 4 T£, 2 tl> 2 WT 2 £, 3 T 3 £, 2 rj> ' 
M-v ~ (A°> WT 2 £, 3 WT 2 WT£, 2 | . (48) 
V T 3 £, 2 iL> WT£, 2 1$ 

The mixing matrix LLv for this mass matrix is given by, 

1 E, 3 

Uv ~ I ^ 1 I • (49) 




Thus we obtain the lepton mixing matrix: 

/ 1 _wi);^\ / 1 0.016 0.014 



u = uliu 




0.021 0.75 -0.66 | . (50) 
6 x TO" 4 0.66 0.75 



Taking (A ) ~ 3 eV, we then obtain suitable masses and mixings for the solution 
of both the solar and atmospheric neutrino problems: 

sin 2 29 atm ~ 1 Ara 2 3 - 1 x 10~ 3 eV 2 sin 2 29 sun - 1 Amf 2 - 6 x 10~ 6 eV 2 (51) 

We did not manage to find an LMA-MSW solution, which is favoured by 
the latest solar neutrino data from Sudbury and SuperKamiokande, using this 
approach. However, during this workshop, Holger, Yasutaka and I constructed 
a promising LMA-MSW solution |9| using the extended version of the AGUT 
model with right-handed neutrinos and the usual see-saw mechanism. 
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1 Introduction 

In this paper, I will update my report 1 1 1 on fermion masses in the Anti-Grand 
Unification Theory (AGUT) at the Bled 2001 workshop. There are two versions of 
the AGUT model based on three family replicated copies of the Standard Model 
(SM) gauge group, selected according to whether or not each family is supple- 
mented by a right-handed neutrino. We note that neither model introduces su- 
persymmetry. In the absence of right-handed neutrinos, the AGUT gauge group 
is Gi = SMG 3 x U(l) f , where SMG = SU(3) x SU(2) x U(1). With the in- 
clusion of three right-handed neutrinos, the AGUT gauge group is extended to 
G2 = (SMG x U(I)b-l) 3 / where the three copies of the SM gauge group are sup- 
plemented by an abelian (B — L) (= baryon number minus lepton number) gauge 
group for each family. In each case, the AGUT gauge group Gi (G2) is the largest 
anomaly free group 1 2 1 transforming the known 45 Weyl f ermions (and the addi- 
tional three right-handed neutrinos for G2) into each other unitarily, which does 
NOT unify the irreducible representations under the SM gauge group. 

Here we present good order of magnitude fits, with four and five adjustable 
parameters respectively, to the quark and lepton masses and mixing angles for 
the two versions of the AGUT model. In each case the fit to the charged fermion 
masses and quark mixings is arranged to essentially reproduce the original three 
parameter AGUT fit Q. It is necessary to introduce a new mass scale into the 
theory, in order to obtain realistic neutrino masses. For the Gi = SMG 3 x U(l)f 
model we introduce a weak isotriplet Higgs field and obtain a vacuum oscillation 
solution to the solar neutrino problem, whereas for the G2 = (SMG x U(I)b-l) 3 
model we introduce the usual see-saw mass scale for the right-handed neutrinos 
and obtain a large mixing angle (LMA) MSW solution. During the last year, fur- 
ther data |3| from the Sudbury Neutrino Observatory have confirmed that the 
LMA-MSW solution is strongly favoured, with the Vacuum Oscillation and LOW 
solutions now allowed at the 3o~ level while the SMA-MSW solution seems to be 
completely ruled out. We refer to 1 4 1 for a recent review of the phenomenology of 
neutrino physics. 
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2 The SMG 3 x U(l) f Model 

The usual SM gauge group is identified as the diagonal subgroup of SMG 3 and 
the AGUT gauge group SMG 3 x U(l)f is broken down to this subgroup by four 
Higgs fields S, W, T and £,. Thus, for example, the SM weak hypercharge y/2 is 
given by the sum of the weak hypercharge quantum numbers \)i/2 for the three 
proto-families: 

2 2 2 2 U 
The spontaneously broken chiral AGUT gauge quantum numbers of the quarks 
and leptons protect the charged fermion masses and generate a mass hierarchy 
for them |5| in terms of the Higgs field vacuum expectation values (VEVs). How- 
ever the VEV of the Higgs field S is taken to be unity in fundamental (Planck) 
mass units. Thus only the VEVs of the other three Higgs fields are used as free 
parameters, in the order of magnitude fit to the effective SM Yukawa coupling 
matrices Yu, Yu arid for the quarks and charged leptons 1 1 [. 

In this model, the large — ~v r mixing required for atmospheric neutrino 
oscillations is generated by introducing a large off-diagonal element, (Ye)23 ~ 
(Ye)33, in the charged lepton Yukawa coupling matrix. This is achieved by intro- 
ducing a new Higgs field 4> with a VEV equal to unity in Planck units and the 
following set of U(l ) gauge charges: 

Q^=3Q t + 2Qw + 4Q T =Q ) -^ ) -y) ( (2) 

Here we express the abelian gauge charges in the model as a charge vector Q = 
( yi /2,u 2 /2,U3/2,Q f ). 

We then have the Yukawa matrices for the quarks: 

/WT 2 £, 2 WT 2 £, W 2 T£,\ / WT 2 f, 2 WT 2 f, T 3 £,\ 

Y u = WT 2 f, 3 WT 2 W 2 T , Y D = WT 2 £, WT 2 T 3 (3) 
V £, 3 1 WT J \ W 3 T W 3 T£,WT/ 

and the charged lepton Yukawa matrix: 

/WT 2 £, 2 W 3 T 2 W 2 T£, 3 \ 
Y E = WT 2 £, 5 WT 2 W£, . (4) 
V W 4 £, 2 W 3 T£, WT / 

We still obtain a good order of magnitude phenomenology for the charged fermion 
masses and quark mixing angles, similar to the original AGUT fit 1 1 1, with the fol- 
lowing VEVs in Planck units: 

<W>= 0.179, <T>= 0.071 <£,>= 0.099. (5) 

The unitary matrix which diagonalises Ye Y^ is then given by 



WE, 4 1 £, 

WE, 5 E, 1 



1 0.05 1.7 x 10~ 4 ^ 
-0.03 0.58 0.81 J (6) 
0.04 -0.81 0.58 
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As we can see from the structure of U E/ we naturally obtain the large u— x mixing 
required for the atmospheric neutrinos. We can now obtain suitable mixing and 
a suitable hierarchy of neutrino masses for a vacuum oscillation solution to the 
solar neutrino problem, by making the following choice of charges for a weak 
iso-triplet Higgs field, A. 



We then have the neutrino mass matrix, 



Mv ~< A > 



WE, 6 WE, 3 WYE 2 ' 
WE, 3 W WYE 
WYE 2 WJE WT 2 £,, 



This has the hierarchy, 



Amf 2 
Am? 3 
Am 2 2 



Ami 3 , 
2J 3 E 3 



7x 10 



-7 



(8) 



(9) 
(10) 



which is just suitable for the atmospheric neutrinos and the vacuum oscillation 
solution to the solar neutrino problem. The electron neutrino mixing is also large 
enough for the vacuum oscillation solution to the solar neutrino problem, as we 
can see from the matrix U^, which diagonalises : 



Uv - 



-M E 3 -Ml 





(ii) 



Hence we have the lepton mixing matrix, 

10 



u = u f E u^ 



WE," 



1 



£.(1- 



l 



\/2(l- 
t(1- 



Al 



T^2(1- 
1 + A 



0.83 2.8 x 10~ 2 0.58 
0.47 0.58 -0.68 
-0.34 0.813 0.49 



(12) 



(13) 



If we take < A >~ 0. 1 8 eV, then we have 



mi - < A > (^-WTE, 2 + ^j^j ~ -1.4 x lO" 5 eV 
m 2 ~ < A > W - 3.: 
m 3 ~ < A > ( WT£, 2 



m 2 ~<A° >W^3.2x 10~ 2 eV 



3 x 10~ 5 eV 



(14) 
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The Am 2 and mixing angle for the solar neutrinos are given by 

Am 2 3 ~7x 1CT 10 eV 2 , sin 2 29 ~ 4U 2 e] U 2 3 -~ 0.93 (15) 

which are compatible with vacuum oscillations for the solar neutrinos. Similarly 
the Am 2 and mixing angle for atmospheric neutrinos are given by, 

Am 2 23 ~~ 1 x 10~ 3 eV 2 , sin 2 2e Qtm - 0.93. (16) 

Hence we can see that we have large (but not maximal) mixing for both the so- 
lar and atmospheric neutrinos. We also note that the CHOOZ electron survival 
probability bound is readily satisfied by U e 2 ~ 0.028 < 0.1 6, which is the relevant 
mixing matrix element since Am 2 2 ~ Am 2 3 ^> Am 2 3 . Thus we obtain a good 
order of magnitude fit (agreeing with the data to within a factor of 2) to the 17 
measured fermion mass and mixing angle variables with just 4 free parameters 
(W, T, E, and A ), but assuming a vacuum oscillation solution to the solar neutrino 
problem. 

We note that we have really only used the abelian gauge quantum num- 
bers to generate a realistic spectrum of fermion masses; the non-abelian repre- 
sentations are determined by imposing the usual SM charge quantisation rule for 
each of the SMG factors in the gauge group. Furthermore two of the U(1)s (more 
precisely two linear combinations of the U(l)s) in the gauge group are sponta- 
neously broken by the Higgs fields S and which have VEVs < S >=< 4> >= 1 . 
Hence these U(1)s play essentially no part in obtaining the spectrum of fermion 
masses and mixings. This means that we can construct a model based on the 
gauge group SMG x U(l )' with the same fermion spectrum as above. However it 
turns out 1 6 [ that some of the quarks and leptons must have extremely large (in- 
teger) U(l)' charges, making this reduced model rather unattractive. Also three 
Higgs fields W, T and £, are responsible for the spontaneous breakdown to the SM 
gauge group, SMG x U(1 )' — > SMG, and they have large relatively prime U(l )' 
charges. So we prefer the better motivated SMG 3 x U(1 )f AGUT model. 

3 The (SMG x U(1 ) B L ) 3 Model 

In this extended AGUT model we introduce a right-handed neutrino and a gauged 
B — L charge for each family with the associated abelian gauge groups U(1) b _l,i 
(i = 1,2, 3). The U(l )f abelian factor of the SMG x U(l )f model in the previous 
section gets absorbed as a linear combination of the B — L charge and the weak 
hypercharge abelian gauge groups for the different families (or generations). It is 
these 6 abelian gauge charges which are responsible for generating the fermion 
mass hierarchy and we list their values in Table^for the 48 Weyl proto-fermions 
in the model. The see-saw scale for the right-handed neutrinos is introduced via 
the VEV of a new Higgs field cf) ss . However, in order to get an LMA-MSW solu- 
tion to the solar neutrino problem, we have to replace 1 8 1 the AGUT Higgs fields 
S and £, by two new Higgs fields p and u>. The abelian gauge quantum numbers 
of the new system of Higgs fields for the (SMG x U(I)b-l) 3 model are given in 
Tabled 
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Table 1. All U ( 1 ) quantum charges for the proto-fermions in the (SMGxII(1) B -l) 3 model. 
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Table 2. All U ( 1 ) quantum charges of the Higgs fields in the (SMG xU(1) B -l) 3 model. 
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As can be seen from Table |5| the fields cu and p have only non- trivial quan- 
tum numbers with respect to the first and second families. This choice of quantum 
numbers makes it possible to express a fermion mass matrix element involving 
the first family in terms of the corresponding element involving the second fam- 
ily, by the inclusion of an appropriate product of powers of p and cu. With the 
system of quantum numbers in Table |3 one can easily evaluate, for a given mass 
matrix element, the numbers of Higgs field VEVs of the different types needed 
to perform the transition between the corresponding left- and right-handed Weyl 
fields. The results of calculating the products of Higgs fields needed, and thereby 
the order of magnitudes of the mass matrix elements in our model, are presented 
in the following mass matrices (where, for clarity, we distinguish between Higgs 
fields and their hermitian conjugates): 
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the up-type quarks: 

((W f ) 




cut)3\vtT 2 tuptW^T 2 cupt(Wt) 2 T^ 
tut) 4 pWtT 2 Wn 2 (Wt) 2 T j (17) 



cu 3 W(Tt) 2 cuptWtTt) 2 tup t T 3N 
cu 2 pW(Tt) 2 W(Tt) 2 T 3 | (18) 
co 2 pW 2 (Tt) 4 W 2 (Tt) 4 VVT 



the charged leptons: 

cu 3 W(Tl) 2 (cut) 3 p 3 W(Tt) 2 (tut) 3 p 3 W 4 (Tt) ; 



(4>v> 




cu 6 (pt) 3 W(Tt) 2 W(Tt) 2 W 4 (Tt) 5 ) (19) 

a,6(pt)3( W t)2 T 4 (wt) 2 T 4 WT 



(cut)3 W t T 2 ( w i)3 p 3\Vij2 (a)t)3 p 3 W 2 (T t) 

(pt)3 W t T 2 w t T 2 W 2 (Tt) 7 | (20) 

( p t)3 (W t)4 T 8 ( W t) 4 T 8 Wf 

and the Majorana (right-handed) neutrinos: 

(pt) 6 T 6 (pt) 3 T 6 (pt) 3 W 3 (Tt)^ 
M R ~(4> SS )| (pt) 3 T 6 T 6 W 3 (Tt) 3 j (21) 

(pt)3 W 3 (T t)3 W 3( T t)3 W 6 (Tt) 12 

Then the light neutrino mass matrix - effective left-left transition Majorana 
mass matrix - can be obtained via the see-saw mechanism 1 7 1 : 

M e ff « Ivt^ Mj, 1 (M°) T . (22) 

with an appropriate renormalisation group running from the Planck scale to the 
see-saw scale and then to the electroweak scale. The experimental quark and lep- 
ton masses and mixing angles in Table |3] can now be fitted, by varying just 5 
Higgs field VEVs and averaging over a set of complex order unity random num- 
bers, which multiply all the independent mass matrix elements. The best fit is 
obtained with the following values for the VEVs: 

(4> ss ) =5.25 x 10 15 GeV , (tu) = 0.244 , (p) = 0.265 , 
(W) =0.157 , (T) =0.0766 , (23) 

where, except for the Higgs field (4> ss ), the VEVs are expressed in Planck units. 
The resulting 5 parameter order of magnitude fit, with an LMA-MSW solution to 
the solar neutrino problem, is presented in Table 

Transforming from tan 2 9 variables to sin 2 20 variables, our predictions for 
the neutrino mixing angles become: 

sin 2 20 Q = 0.66 , sin 2 20 atm = 0.96 , sin 2 20 chooz = 0.1 1 . (24) 
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Table 3. Best fit to conventional experimental data. All masses are running masses at 1 GeV 
except the top quark mass which is the pole mass. 





Fitted 
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4.4 MeV 


4 MeV 


TUd 


4.3 MeV 


7 MeV 


m e 


1 .6 MeV 


0.5 MeV 


m c 


0.64 GeV 


1 .4 GeV 


m s 


295 MeV 


T /\/\ A iff T7 

200 MeV 




111 MeV 


1 /\ r" A iff T 7 

1 05 MeV 


M, 


202 GeV 


1 80 GeV 


m.b 


5.7 GeV 


6.3 GeV 




1 .46 GeV 


1 .78 GeV 


Vus 


0.11 


0.22 


V cb 


0.026 


0.041 


V ub 


0.0027 


0.0035 


Am Q 


9.0 x lO - * eV 2 


5.0 x )0~ b eV 2 


Am 2 lm 


1.7 x 10~ 3 eV 2 


2.5 x 10~ 3 eV 2 


tan 2 e 


0.26 


0.34 


tan 2 9 a tm 


0.65 


1.0 


tan 2 6 c hooz 


2.9 x 10~ 2 


$2.6 x 10~ 2 



Note that our fit to the CHOOZ mixing angle lies close to the 2o~ Confidence Level 
experimental bound. We also give here our predicted hierarchical left-handed 
neutrino masses (mi) and the right-handed neutrino masses (Mi) with mass eigen- 
state indices (i = 1 , 2, 3): 

mi = 1.4 x 1(T 3 eV , M, = 1 .0 x 10 6 GeV , (25) 

m 2 = 9.6 x ICr 3 eV , M 2 = 6.1 x 10 9 GeV , (26) 
m 3 = 4.2 x 1CT 2 eV , M 3 = 7.8 x 10 9 GeV . (27) 
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